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Abstract. Wg provide an introduction to selected recent advances in the 
mathematical understanding of Einstein's theory of gravitation. 
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1. Introduction 

Mathematical general relativity is, by now, a well-established vibrant branch of 
mathematics. It ties fundamental problems of gravitational physics with beautiful 
questions in mathematics. The object is the study of manifolds equipped with a 
Lorentzian metric satisfying the Einstein field equations. Some highlights of its his- 
tory include the discovery by Choquet-Bruhat of a well posed Cauchy problem [178], 
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subsequently globalized by Choquet-Bruliat and Geroch [93], the singularity the- 
orems of Penrose and Hawking [201,317], the proof of the positive mass theorem 
by Schocn and Yau [344], and the proof of stability of Minkowski space-time by 
Christodoulou and Klainerman [107]. 

There has recently been spectacular progress in the field on many fronts, includ- 
ing the Cauchy problem, stability, cosmic censorship, construction of initial data, 
and asymptotic behaviour, many of which will be described here. Mutual bene- 
fits are drawn, and progress is being made, from the interaction between general 
relativity and geometric analysis and the theory of elliptic and hyperbolic partial 
differential equations. The Einstein equation shares issues of convergence, collapse 
and stability with other important geometric PDEs, such as the Ricci flow and 
the mean curvature flow. Steadily growing overlap between the relevant scientific 
communities can be seen. For all these reasons it appeared timely to provide a 
mathematically oriented reader with an introductory survey of the field. This is 
the purpose of the current work. 

In Section 2 we survey the Lorentzian causality theory, the basic language for de- 
scribing the structure of space-times. In Section 3 the reader is introduced to black 
holes, perhaps the most fascinating prediction of Einstein's theory of gravitation, 
and the source of many deep (solved or unsolved) mathematical problems. In Sec- 
tion 4 the Cauchy problem for the Einstein equations is considered, laying down the 
foundations for a systematic construction of general space-times. Section 5 exam- 
ines initial data sets, as needed for the Cauchy problem, and their global properties. 
In Section 6 we discuss the dynamics of the Einstein equations, including questions 
of stability and predictability; the latter question known under the baroque name 
of "strong cosmic censorship" . Section 7 deals with trapped and marginally trapped 
surfaces, which signal the presence of black holes, and have tantalizing connections 
with classical minimal surface theory. The paper is sprinkled with open problems, 
which are collected in Appendix A. 

2. Elements of Lorentzian geometry and causal theory 

2.1. Lorentzian manifolds. In general relativity, and related theories, the space 
of physical events is represented by a Lorentzian manifold. A Lorentzian manifold 
is a smooth (Hausdorff, paracompact) manifold ^ ~ of dimension n + 1, 

equipped with a Lorentzian metric g. A Lorentzian metric is a smooth assignment 
to each point p e ^ of a symmetric, nondegenerate bilinear form on the tangent 

space Tj^^ of signature ( h •••+). Hence, if {eo, ei, e„} is an orthonormal 

basis for with respect to g, then, perhaps after reordering the basis, the matrix 
[g{ei, Ej)] equals diag (—1, -1-1, -1-1). A vector v = J2 v^ea then has 'square norm', 



which can be positive, negative or zero. This leads to the causal character of vectors, 
and indeed to the causal theory of Lorentzian manifolds, which we shall discuss in 
Section 2.3. 

On a coordinate neighborhood {U, x°')= {U, ) the metric g is com- 

pletely determined by its metric component functions on U , gap '■= 5(af^, af?)' 
< Q,/3 < n: For v = V^^^w = w^^ ^ € Tp^, P € U, g{v,w) = gapv^w^. 
(Here we have used the Einstein summation convention: If, in a coordinate chart, 
an index appears repeated, once up and once down, then summation over that 



(2.1) 
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index is implied.) Classically the metric in coordinates is displayed via the "line 
element", ds"^ = gapdx^dx^ . 

The prototype Lorcntzian manifold is Minkowski space R^'", the space-time of 
special relativity. This is equipped with the Minkowski metric, which, with 

respect to Cartesian coordinates {x'^,x^, ...,x"), is given by 

ds^ = -{dx°f + [dx^Y + ■■■ + {dx"f . 

Each tangent space of a Lorentzian manifold is isometric to Minkowski space, and 
in this way the local accuracy of special relativity is built into general relativity. 

Every Lorentzian manifold (or, more generally, pseudo-Riemannian manifold) 
(^""'"^,(7) comes equipped with a Levi-Civita connection (or covariant differentia- 
tion operator) V that enables one to compute the directional derivative of vector 
fields. Hence, for smooth vector fields X,Y G 3i{^), VxY G X{^) denotes the co- 
variant derivative of Y in the direction X. The Levi-Civita connection is the unique 
connection V on (^"+^,3) that is (i) symmetric (or torsion free), i.e., that satisfies 
VxY - VyX = [X,Y] for aU X,Y e X(M), and (ii) compatible with the metric, 
i.e. that obeys the metric product rule, X{g{Y, Z)) = g{WxY, Z) + g{Y, WxZ), for 
all X,Y,Z €X{M). 

In a coordinate chart (U, x") , one has, 

(2.2) VxY = {xiY^ + r^^x"y^)a^ , 

where X", Y" are the components of X and Y, respectively, with respect to the co- 
ordinate basis da = g§^, and where the Ejj^^'s are the classical Christoffcl symbols, 
given in terms of the metric components by, 

(2.3) r^^ = Ig^'idpga. + d^gpu - d^g^p) . 

Note that the coordinate expression (2.2) can also be written as, 

(2.4) VxY = X^VaY^^d^, , 
where V^Y^ (often written classically as Y^.a) is given by, 

(2.5) v„y^ = a^y^ r^^F''. 

We shall feel free to interchange between coordinate and coordinate free no- 
tations. The Levi-Civita connection V extends in a natural way to a covariant 
differentiation operator on all tensor fields. 

The Riemann curvature tensor of (^"+^,5) is the map R : X{M) x X{M) x 
X{M) X(M), {X, Y, Z) R{X, Y)Z, given by 

(2.6) R{X, Y)Z = Vx^yZ - Vy'^xZ - V[x,y]^ . 

This expression is linear in X, y, Z € X(^) with respect to C°°(>^). This implies 
that R is indeed tensorial, i.e., that the value of RiX, Y)Z at p £ M depends only 
on the value of X, Y, Z at p. 

Equation (2.6) shows that the Riemann curvature tensor measures the extent 
to which covariant differentiation fails to commute. This failure to commute may 
be seen as an obstruction to the existence of parallel vector fields. By Riemann's 
theorem, a Lorentzian manifold is locally Minkowskian if and only if the Riemann 
curvature tensor vanishes. 

The components R^-^ap of the Riemann curvature tensor i? in a coordinate chart 
(f7, a;")are determined by the equations, R{dcndp)d^ — R^-^apd^. Equations (2.2) 
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and (2.6) then yield the following explicit formula for the curvature components in 
terms of the Christoffel symbols, 

(2.7) = 9„r!;^ - dpT^^^ + r!;^r(;„ - v^j^^^ . 

The Ricci tensor, Ric, is a bilinear form obtained by contraction of the Riemann 
curvature tensor, i.e., its components = 'Ric{d^,du) are determined by tracing, 
R^J.u = R" nay Symmetries of the Riemann curvature tensor imply that the Ricci 
tensor is symmetric, R^u = Rly^l■ By tracing the Ricci tensor, wc obtain the scalar 
curvature, R = g^'^ R^u, where g^^'^ denotes the matrix inverse to g^^- 

2.2. Einstein equations. The Einstein equation (with cosmological constant A), 
the field equation of general relativity, is the tensor equation, 

(2.8) Ric-ii?g + A5 = 87rr, 

where T is the energy-momentum tensor. (See, e.g.. Section 6.5.2 for an example of 
an energy-momentum tensor.) When expressed in terms of coordinates, the Einstein 
equation becomes a system of second order equations for the metric components g^i/ 
and the nongravitational field variables introduced through the energy-momentum 
tensor. We say that space-time obeys the vacuum Einstein equation if it obeys the 
Einstein equation with T = 0. 

The Riemann curvature tensor has a number of symmetry properties, one of 
which is the so-called first Bianchi identity: 

RaP^S + Ra-ySp + RaSp-y = . 

The curvature tensor also obeys a differential identity known as the second Bianchi 
identity: 

(2.9) ^(jRafSjS + '^aRpa',S + '^pRaa',S = ■ 

When twice contracted, (2.9) yields the following divergence identity: 

(2.10) v„ (^i?"^ - l^"^^ = . 

This plays a fundamental role in general relativity, as, in particular, it implies, in 
conjunction with the Einstein equation, local conservation of energy, V^T"^ = 0. It 
also plays an important role in the mathematical analysis of the Einstein equations; 
see Section 4 for further discussion. 

2.3. Elements of causal theory. Many concepts and results in general relativity 
make use of the causal theory of Lorentzian manifolds. The starting point for causal 
theory is the causal classification of tangent vectors. Let (^"+^, g) be a Lorentzian 
manifold. A vector v E is timelike (resp., spacelikc, null) provided g(y, u) < 
(resp., g{v,v) > 0, g{v,v) = 0). The collection of null vectors forms a double cone 
Vp in (recall (2.1)), called the null cone at p; see Figure 2.1. 

The timelike vectors at p point inside the null cone and the spacelike vectors 
point outside. We say that v g Tp^ is causal if it is timelike or null. We define the 
length of causal vectors as |u| = \/—g{v, v). Causal vectors v,w G Tp^ that point 
into the same half-cone of the null cone Vp obey the reverse triangle inequality, 
b + > \v\ -f \w\. Geometrically, this is the source of the twin paradox. 

These notions of causality extend to curves. Let 7 : / — t — 7(t), be a 
smooth curve in then 7 is said to be timelike (resp., spacelike, null, causal) 
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past pointing timelike 



past pointing null 





future pointing null 



future pointing timelike 



Figure 2.1. The light cone at p. 



provided each of its velocity vectors j'{t) is timelike (resp., spacelike, null, causal). 
Heuristically, in accordance with relativity, information flows along causal curves, 
and so such curves are the focus of attention in causal theory. The notion of a 
causal curve extends in a natural way to piecewise smooth curves, and we will 
normally work within this class. As usual, we define a geodesic to be a curve 
t — ^ j(t) of zero covariant acceleration, = 0. Since geodesies 7 are constant 

speed curves (5(7', 7') = const.), each geodesic in a Lorentzian manifold is either 
timelike, spacelike or null. 

The length of a causal curve 7 : [a, b] — s- is defined as 



If 7 is timelike one can introduce an arc length parameter along 7. In general 
relativity, a timelike curve corresponds to the history of an observer, and arc length 
parameter, called proper time, corresponds to time kept by the observer. Using 
the existence and properties of geodesically convex neighborhoods [307] one can 
show that causal geodesies arc locally maximal (i.e., locally longest among causal 
curves). 

Each null cone Vp consists of two half-cones, one of which may designated as 
the future cone, and the other as the past cone at p. If the assignment of a past 
and future cone at each point of ^ can be carried out in a continuous manner 
over ^ then ^ is said to be time-orientable. There are various ways to make 
the phrase "continuous assignment" precise, but they all result in the following 
fact: A Lorentzian manifold (^"'*'^, 5) is time-orientable if and only if it admits 
a smooth timelike vector field Z. If ^ is time-orientable, the choice of a smooth 
timelike vector field Z fixes a time orientation on For any p g a causal 
vector V G is future directed (resp. past directed) provided g{v, Z) < (resp. 

g{v,Z) > 0). Thus, v is future directed if it points into the same null half cone 
at p as Z . We note that a Lorentzian manifold that is not time-orientable always 
admits a double cover that is. By a space-time we mean a connected time-oriented 
Lorentzian manifold (./#"+^,(7). Henceforth, we restrict attention to space-times. 

2.3.1. Past and futures. Let (^, g) be a space-time. A timelike (resp. causal) 
curve 7 : / — s- ^ is said to be future directed provided each tangent vector ^'{t), 



L{j) = Length of 7 
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t G /, is future directed. (Past- directed timelike and causal curves are defined in a 
time-dual manner.) I^{p), the timelike future of p S is the set consisting of 
all points q € ^ for which there exists a future directed timelike curve from p to 
q. J'^{p), the causal future of p S is the set consisting oi p and all points q for 
which there exists a future directed causal curve from p to q. In Minkowski space 
R-*^'" these sets have a simple structure: For each p € M^'", dl'^{p) = J'^{p) \ I~^{p) 
is the future cone at p generated by the future directed null rays emanating from 
p. I^{p) consists of the points inside the cone, and J^{p) consists of the points 
on and inside the cone. In general, the curvature and topology of space-time can 
strongly influence the structure of these sets. 

Since a timelike curve remains timelike under small smooth perturbations, it is 
heuristically clear that the sets {p) are in general open; a careful proof makes use 
of properties of geodesically convex sets. On the other hand the sets J'^ip) need 
not be closed in general, as can be seen by considering the space-time obtained by 
removing a point from Minkowski space. 

It follows from variational arguments that, for example, if q £ I^{p) a-nd r G 
J^{q) then r G I^{p)- This and related claims are in fact a consequence of the 
following fundamental causality result [307]. 

Proposition 2.1. If q G J+(p)\/+(p), i.e., if q is in the causal future of p but not 
in the timelike future of p then any future directed causal curve from p to q must 
be a null geodesic. 

Given a subset S C I^{S), the timelike future of S, consists of all points 
q G ^ for which there exists a future directed timelike curve from a point in 5* to 
q. J~^{S), the causal future of S consists of the points of S and all points q G ./# 
for which there exists a future directed causal curve from a point in S to q. Note 
that I'^iS) = Upes ^'^(p)- Hence, as a union of open sets, I~^{S) is always open. 

The timelike and causal pasts I~{p), J~{p), I~{S), J~{S) are defined in a time 
dual manner in terms of past directed timelike and causal curves. It is sometimes 
convenient to consider pasts and futures within some open subset U of For 
example, I~^{p,U) denotes the set consisting of all points g G C/ for which there 
exists a future directed timelike from p to q contained in U. 

Sets of the form dI^{S) are called achronal boundaries^ and have nice structural 
properties: They are achronal Lipschitz hypersurfaces, ruled, in a certain sense, by 
null geodesies [307]. (A set is achronal if no two of its points can be joined by a 
timelike curve.) 

2.3.2. Causality conditions. A number of results in Lorentzian geometry and gen- 
eral relativity require some sort of causality condition. It is perhaps natural on 
physical grounds to rule out the occurrence of closed timelike curves. Physically, 
the existence of such a curve signifies the existence of an observer who is able 
to travel into his/her own past, which leads to variety of paradoxical situations. 
A space-time ^ satisfies the chronology condition provided there are no closed 
timelike curves in It can be shown that all compact space-times violate the 
chronology condition, and for this reason compact space-times have been of limited 
interest in general relativity. 

A somewhat stronger condition than the chronology condition is the causality 
condition. A space-time ^ satisfies the causality condition provided there are 
no closed (nontrivial) causal curves in A slight weakness of this condition is 
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that there are space-times which satisfy the causahty condition, but contain causal 
curves that are "ahnost closed", see e.g. [204, p. 193]. 

It is useful to have a condition that rules out "almost closed" causal curves. A 
space-time is said to be strongly causal at p € y/^ provided there arc arbitrarily 
small neighborhoods U oip such that any causal curve 7 which starts in, and leaves, 
U never returns to U. ^ is strongly causal if it is strongly causal at each of its 
points. Thus, heuristically speaking, ^ is strongly causal provided there are no 
closed or "almost closed" causal curves in Strong causality is the "standard" 
causality condition of space-time geometry, and although there are even stronger 
causality conditions, it is sufficient for most applications. A very useful fact about 
strongly causal space-times is the following: If ^ is strongly causal then any 
future (or past) inextendible causal curve 7 cannot be "imprisoned" or "partially 
imprisoned" in a compact set. That is to say, if 7 starts in a compact set K, it 
must eventually leave K for good. 

We now come to a fundamental condition in space-time geometry, that of global 
hyperbolicity. Mathematically, global hypcrbolicity is a basic 'niceness' condition 
that often plays a role analogous to geodesic completeness in Ricmannian geom- 
etry. Physically, global hypcrbolicity is connected to the notion of strong cosmic 
censorship, the conjecture that, generically, space-time solutions to the Einstein 
equations do not admit naked (i.e., observable) singularities; see Section 6.2 for 
further discussion. 

A space-time ^ is said to be globally hyperbolic provided: 

(1) ^ is strongly causal. 

(2) (Internal Compactness) The sets J+(p)n {q) are compact for allp, q G 

Condition (2) says roughly that has no holes or gaps. For example Minkowski 
space M^'" is globally hyperbolic but the space-time obtained by removing one 
point from it is not. Leray [259] was the first to introduce the notion of global 
hypcrbolicity (in a somewhat different, but equivalent form) in connection with his 
study of the Cauchy problem for hyperbolic PDEs. 

We mention a couple of basic consequences of global hypcrbolicity. Firstly, glob- 
ally hyperbolic space-times are causally simple, by which is meant that the sets 
J^{A) are closed for all compact A C This fact and internal compactness 

implies that the sets J^(^) n J^{B) are compact, for all compact A,Bc ^ . 

Analogously to the case of Riemannian geometry, one can learn much about the 
global structure of space-time by studying its causal geodesies. Global hyperbolicity 
is the standard condition in Lorentzian geometry that guarantees the existence of 
maximal timelike geodesic segments joining timelike related points. More precisely, 
one has the following. 

Proposition 2.2. // is globally hyperbolic and q G /^(p), then there exists 
a maximal timelike geodesic segment 7 from p to q ( where by maximal, we mean 
Lij) > L((7) for all future directed causal curves a from p to q). 

Contrary to the situation in Ricmannian geometry, geodesic completeness does 
not guarantee the existence of maximal segments, as is well illustrated by anti-de 
Sitter space, see e.g. [31]. 

Global hypcrbolicity is closely related to the existence of certain 'ideal initial 
value hypersurfaces', called Cauchy (hyper) surf aces. There are slight variations in 
the literature in the definition of a Cauchy surface. Here we adopt the following 
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definition: A Caucliy surface for a space-time ^ is a subset S that is met exactly 
once by every inextendible causal curve in ^ . It can be shown that a Cauchy 
surface for ^ is necessarily a C° (in fact, Lipschitz) hypersurface in . Note also 
that a Cauchy surface is acausal, that is, no two of its points can be joined by a 
causal curve. The following result is fundamental. 

Proposition 2.3 (Geroch [192]). ^ is globally hyperbolic if and only if admits 
a Cauchy surface. If S is a Cauchy surface for then ^ is homeomorphic to 
K X S". 

With regard to the implication that global hyperbolicity implies the existence of 
a Cauchy surface, Geroch, in fact, proved something substantially stronger. (We 
will make some comments about the converse in Section 2.3.3.) A time function 
on ^ is a C° function t on ^ such that t is strictly increasing along every future 
directed causal curve. Geroch established the existence of a time function t all of 
whose level sets t = fo, G M, are Cauchy surfaces. This result can be strengthened 
to the smooth category. By a smooth time function we mean a smooth function 
t with everywhere past pointing timelike gradient. This implies that t is strictly 
increasing along all future directed causal curves, and that its level sets are smooth 
spacelike^ hypersurfaccs. It has been shown that a globally hyperbolic space-time 
admits a smooth time function all of whose levels sets are Cauchy surfaces [46,352]. 
In fact, one obtains a diffeomorphism ^ w M x S", where the M-factor corresponds 
to a smooth time function, such that each slice St = {t} x S", t G K, is a Cauchy 
surface. 

Given a Cauchy surface S to begin with, to simply show that ^ is homeomorphic 
to M x iS, consider a complete timelike vector field Z on ^ and observe that each 
integral curve of Z, when maximally extended, meets 5 in a unique point. This 
leads to the desired homeomorphism. (If S is smooth this will be a diffeomorphism.) 
In a similar vein, one can show that any two Cauchy surfaces are homeomorphic. 
Thus, the topology of a globally hyperbolic space-time is completely determined by 
the common topology of its Cauchy surfaces. 

The following result is often useful. 

Proposition 2.4. Let ^ be a space-time. 

(1) If S is a compact acausal C'^ hypersurface and ^ is globally hyperbolic then 
S must be a Cauchy surface for ^ . 

(2) If t is a smooth time function on ^ all of whose level sets are compact, 
then each level set is a Cauchy surface for ^ , and hence ^ is globally 
hyperbolic. 

We will comment on the proof shortly, after Proposition 2.6. 

2.3.3. Domains of dependence. The future domain of dependence of an acausal set 
S is the set &^{S) consisting of all points p £ ^ such that every past inextendible 
causal curve^ from p meets S . In physical terms, since information travels along 
causal curves, a point in ^+(5") only receives information from S. Thus, in princi- 
ple, Sl^ {S) represents the region of space-time to the future of 5* that is predictable 



A hypersurface is called spacelike if the induced metric is Ricmannian; see Section 2.4. 
^We note that some authors use past inextendible timelike curves to define the future domain 
of dependence, which results in some small differences in certain results. 
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Figure 2.2. Examples of domains of dependence and Cauchy 
horizons. 

from S. J^^{S), the future Cauchy horizon of S, is defined to be the future bound- 
ary of ^+(5); in preeise terms, J^f+{S) = {p & &+{S) : I+ip) f] &+{S) = 0}. 
Physically, ^^+(5*) is the future limit of the region of space-time predictable from 
S . Some examples of domains of dependence, and Cauchy horizons, can be found 
in Figure 2.2. 

It follows almost immediately from the definition that is achronal. In 

fact, Cauchy horizons have structural properties similar to achronal boundaries, as 
indicated in the following. 

Proposition 2.5. Let S he an acausal subset of a space-time ^ . Then (S)\S , 
if nonempty, is an achronal C*^ hypersurface of ^ ruled by null geodesies, called 
generators, each of which either is past inextendible in ^ or has past end point on 
S. 

The past domain of dependence i^~(S') of S, and the past Cauchy horizon 
{S) of 5*, are defined in a time-dual manner. The total domain of depen- 
dence &{S) and the total Cauchy horizon are defined respectively as, 
S){S) = S>+{S) U ^-{S) and M'iS) = M'+{S) U J^-{S) 

Domains of dependence may be used to characterize Cauchy surfaces. In fact, 
it follows easily from the definitions that an acausal subset S C ^ is a Cauchy 
surface for ^ if and only if &{S) = ^. Using the fact that dSi{S) ^ (5), we 
obtain the following. 

Proposition 2.6. Let S be an acausal subset of a space-time ./#. Then, S is a 
Cauchy surface for if and only if S){S) = if and only if ,yif{S) ~ 0. 

Part 1 of Proposition 2.4 can now be readily proved by showing, with the aid of 
Proposition 2.5, that J^{S) = 0. Indeed if Jff~^{S) ^ then there exists a past 
inextendible null geodesic r\ C with future end point p imprisoned in the 

compact set (S) H J" (p) which, as already mentioned, is not possible in strongly 
causal space-times. Part 2 is proved similarly; compare [69,187]. 

The following basic result ties domains of dependence to global hyperbolicity. 

Proposition 2.7. Let S C ^ be acausal. 

(1) Strong causality holds at each point of int !^(S). 
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(2) Internal compactness holds on int 2i{S), i.e., for allp, q G int 2i{S), J+(p)n 
(q) is compact. 

Propositions 2.6 and 2.7 immediately imply that if S* is a Cauchy surface for a 
space-time ^ then ^ is globally hyperbolic, as claimed in Proposition 2.3. 

2.4. Submanifolds. In addition to curves, one may also speak of the causal char- 
acter of higher dimensional submanifolds. Let be a smooth submanifold of a 
space-time {^,g). For p & V, we say that the tangent space TpV is spacelike 
(resp. timelike, null) provided g restricted to TpV is positive definite (resp., has 
Lorentzian signature, is degenerate). Then V is said to be spacelike (resp., time- 
like, null) provided each of its tangent spaces is spacelike (resp., timelike, null). 
Hence if V is spacelike (resp., timelike) then, with respect to its induced metric, 
i.e., the metric g restricted to the tangent spaces of V, F is a Riemannian (resp., 
Lorentzian) manifold. 

3. Stationary black holes 

Perhaps the first thing which comes to mind when general relativity is men- 
tioned arc black holes. These are among the most fascinating objects predicted by 
Einstein's theory of gravitation. Although they have been studied for years, ^ they 
still attract tremendous attention in the physics and astrophysics literature. It is 
seldom realized that, in addition to Einstein's gravity, several other field theories 
are known to possess solutions which exhibit black hole properties, amongst which:^ 

• The "dumb holes" , arising in Euler equations, which are the sonic counter- 
parts of black holes, first discussed by Unruh [367] (compare [111]). 

• The "optical" ones - the black-hole-type solutions arising in the theory of 
moving dielectric media, or in non-linear electrodynamics [258,304]. 

The numerical study of black holes has become a science in itself, see [77, 322] and 
references therein. The evidence for the existence of black holes in our universe is 
growing [172,194,239,289,381]. Reviews of, and further references to, the quantum 
aspects of black holes can be found in [18,68,171,221,312,370]. 

In this section we focus attention on stationary black holes that are solutions 
of the vacuum Einstein equations with vanishing cosmological constant, with one 
exception: the static electro- vacuum Majumdar-Papapetrou solutions, an example 
of physically significant multiple black holes. By definition, a stationary space-time 
is an asymptotically flat space-time which is invariant under an action of R by 
isometrics, such that the associated generator — referred to as Killing vector — 
is timelike^ in the asymptotically flat region. These model steady state solutions. 
Stationary black holes are the simplest to describe, and most mathematical results 
on black holes, such as the uniqueness theorems discussed in Section 3.9, concern 
those. It should, however, be kept in mind that one of the major open problems 

■^The reader is referred to the introduction to [83] for an excellent concise review of the history 
of the concept of a black hole, and to [82,237] for a more detailed one. 

^An even longer list of models and submodels can bo found in [20], see also [19,305]. 

^In fact, in the literature it is always implicitly assumed that the stationary Killing vector K 
is uniformly timelike in the asymptotic region M^xt', by this we mean that g{K,K) < — e < 
for some e and for all r large enough. This uniformity condition excludes the possibility of a 
timelike vector which asymptotes to a null one. This involves no loss of generality in well behaved 
space-times: indeed, uniformity always holds for Killing vectors which are timelike for all large 
distances if the conditions of the positive energy theorem are met [34,138]. 
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ill mathematical relativity is the understanding of the dynamical behavior of black 
hole space-times, about which not much is yet known (compare Section 6.6). 

3.1. The Schwarzschild metric. The simplest stationary solutions describing 
compact isolated objects are the spherically symmetric ones. According to BirkhofF's 
theorem [49], any {n+ l)-dimensional, n > 3, spherically symmetric solution of the 
vacuum Einstein equations belongs to the family of Schwarzschild metrics, param- 
eterized by a mass parameter to: 

(3.1) g ^ ~V^dt^ + V~^dr^ +r^dn^ , 

(3.2) = 1- ^ , t e M , r e (2m, oo) . 

Here dOF' denotes the metric of the standard (n — l)-sphcrc. (This is true without 
assuming stationarity.) 

From now on we assume n = 3, though identical results hold in higher dimen- 
sions. 

We will assume 

TO > , 

because m < leads to metrics which are called "nakedly singular" ; this deserves 
a comment. For Schwarzschild metrics we have 

(3.3) Ro.p.sR'-'^' ^^""^ 



in dimension 3 -f 1, which shows that the geometry becomes singular as r = is 
approached; this remains true in higher dimensions. As we shall see shortly, for 
TO > the singularity is "hidden" behind an event horizon, while this is not the 
case for TO < 0. 

One of the first features one notices is that the metric (3.1) is singular as r = 2to 
is approached. It turns out that this singularity is related to an unfortunate choice 
of coordinates (one talks about "a coordinate singularity"); the simplest way to see 
this is to replace t hy a new coordinate v defined as 

(3.4) v = t + f{r), r = ^ 



leading to 



This brings g to the form 

(3.5) g = -{l~—)dv^ + 2dvdr + r^dn^ 



V 2TOln(r — 2to) . 

2to, 



r 



Wc have det g = sin with all coefficients of g smooth, which shows that g is 
a well-defined Lorentzian metric on the set 

(3.6) weR, re(0,oo). 

More precisely, (3.5)-(3.6) provides an analytic extension of the original space-time 
(3.1). 

We could have started immediately from the form (3.5) of 5, which would have 
avoided the lengthy discussion of the coordinate transformation (3.4). However, the 
form (3.1) is the more standard one. Furthermore, it makes the metric g manifestly 
asymptotically flat (see Section 3.6); this is somewhat less obvious to an untrained 
eye in (3.5). 
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It is easily seen that the region {r < 2m} for the metric (3.5) is a black hole 
region, in the sense that 

(3.7) observers, or signals, can enter this region, but can never leave it. 

In order to see that, recall that observers in general relativity always move on 
future directed timelike curves, that is, curves with timelike future directed tan- 
gent vector. For signals, the curves are causal future directed. Let, then, 7(5) = 
(w(s), r{s), 0(s), (p{s)) be such a timelike curve; for the metric (3.5) the timelikeness 
condition (7(7,7) < reads 

2m 

-(1 )i)2 + 2vf + r'^(f + sin^ d^'^) < . 

r 

This implies 

/ 2777 \ 

vy - (1 -jv + 2fj < . 

It follows that V does not change sign on a timelike curve. The usual choice of time 
orientation corresponds to ii > on future directed curves, leading to 

,^ 2771, . 

-(1 )v + 2f <0 . 

r 

For r < 2777 the first term is non-negative, which enforces /• < on all future 
directed timelike curves in that region. Thus, r is a strictly decreasing function 
along such curves, which implies that future directed timelike curves can cross the 
hypersurface {r = 2m} only if coming from the region {r > 2777}. This motivates 
the name black hole event horizon for {r = 2m, v G M}. The same conclusion (3.7) 
applies for causal curves: it suffices to approximate a causal curve by a sequence of 
timelike ones. 

The transition from (3.1) to (3.5) is not the end of the story, as further extensions 
are possible. For the metric (3.1) a maximal analytic extension has been found 
independently by Kruskal [254], Szekeres [360], Synge [358] and Fronsdal [186]; for 
some obscure reason neither Synge nor Fronsdal are almost never mentioned in this 
context. This extension is depicted^ in Figure 3.1. The region / there corresponds 
to the space-time (3.1), while the extension just constructed corresponds to the 
regions / and //. 

The Kruskal-Szekeres extension is singled out by being maximal in the class 
of vacuum, analytic, simply connected space-times, with all maximally extended 
geodesies 7 either complete, or with the curvature scalar Rais-ysR"^^^ diverging 
along 7 in finite affine time. 

An alternative convenient representation of the Schwarzschild metrics, which 
makes the space-part of g manifestly conformally flat, is given by 



are grateful to J. -P. Nicolas for allowing us to use his figure from [302]. 
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Figure 3.1. The Cartcr-Pcnrose diagram^ for the Kruskal- 
Szekeres space-time with mass M . There are actuaUy two asymp- 
totically flat regions, with corresponding event horizons defined 
with respect to the second region. Each point in this diagram rep- 
resents a two-dimensional sphere, and coordinates are chosen so 
that light-cones have slopes plus minus one. 



3.2. Rotating black holes. Rotating generalizations of the Schwarzschild metrics 
arc given by the family of Kerr metrics^ parameterized by a mass parameter m and 
an angular momentum parameter a. One explicit coordinate representation of the 
Kerr metric is 

(3.9) g = -{l- j dv"^ + 2drdv + Y.dB'^ - 2a sin^ 

+ - sm^ Odcj)^ d(j)dv , 



where 



S = + cos^ 



Note that (3.9) reduces to the Schwarzschild solution in the representation (3.5) 
when a = 0. The reader is referred to [80,308] for a thorough analysis. All Kerr 
metrics satisfying 

2 \ 2 
m > a 

provide, when appropriately extended, vacuum space-times containing a rotating 
black hole. Higher dimensional analogues of the Kerr metrics have been constructed 
by Myers and Perry [297]. 

A fascinating class of black hole solutions of the 4 + 1 dimensional station- 
ary vacuum Einstein equations has been found by Emparan and Rcall [169] (see 
also [122,123,168,170,171,321]). The solutions, called hlack rings, are asymptot- 
ically Minkowskian in spacelike directions, with an event horizon having S*^ x 5^ 
cross-sections. The "ring" terminology refers to the factor in 5^ x 5^. 
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3.3. Killing horizons. Before continuing some general notions are in order. By 
definition, a Killing field is a vector field the local flow of which preserves the metric. 
Killing vectors are solutions of the over-determined system of Killing equations 

(3.10) Vc^Xp + V^Xa = . 

One of the features of the metric (3.1) is its stationarity, with Killing vector field 
X ^ dt'. As already pointed out, a space-time is called stationary if there ex- 
ists a Killing vector field X which approaches dt in the asymptotically flat region 
(where r goes to oo, see Section 3.6 for precise definitions) and generates a one 
parameter group of isometrics. A space-time is called static if it is stationary and 
if the distribution of hypcrplanes orthogonal to the stationary Killing vector X is 
integrable. 

A space-time is called axisymmetric if there exists a Killing vector field Y which 
generates a one parameter group of isometrics and which behaves like a rotation: 
this property is captured by requiring that all orbits be 27r-periodic, and that the 
set {Y = 0}, called the axis of rotation, be non-empty. 

Let X be a Killing vector field on {^,g), and suppose that ^ contains a 
null hypersurface (see Sections 2.4 and 7.1) ^ = ,Ao(X) which coincides with a 
connected component of the set 

.yViX) := {pe^\ g{Xp, Xp) = , ^ 0} , 

with X tangent to -y^Q. Then is called a Killing horizon associated to the Killing 
vector X. The simplest example is provided by the "boost Killing vector field" 

(3.11) K^zdt+td, 

in four-dimensional Minkowski space-time M}''^: ,jV{K) has four connected compo- 
nents 

J^{K),s:={t = ez ,5t>{)} , e,(5e{±l}. 

The closure JV{K) of .vV{K) is the set {\t\ = |z|}, which is not a manifold, because 
of the crossing of the null hyperplanes {t = ±z} at t = 2 = 0. Horizons of this type 
are referred to as bifurcate Killing horizons. 

A very similar behavior is met in the extended Schwarzschild space-time: the set 
{r = 2m} is a null hypersurface the Schwarzschild event horizon. The stationary 
Killing vector X = dt extends to a Killing vector X which becomes tangent to and 
null on S' in the extended space-time, except at the "bifurcation sphere" right in 
the middle of Figure 3.1, where X vanishes. 

A last noteworthy example in Minkowski space-time M}'^ is provided by the 
Killing vector 

(3.12) X = ydt + tdy + xdy - yd^ = ydt + [t + x)dy - yd^ ■ 

Thus, X is the sum of a boost ydt + tdy and a rotation xdy — ydx- Note that X 
vanishes if and only if 

y ~ t + X = , 

which is a two-dimensional null submanifold of M.^'^. The vanishing set of the 
Lorcntzian length of X , 

giX,X) = (f + a;)2 = , 
is a null hyperplane in R^-^. It follows that, e.g., the set 

{t + x^O , ?/>0,t>0} 
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is a Killing horizon with respect to two different Killing vectors, the boost Killing 
vector xdt + tdx, and the Killing vector (3.12). 

3.3.1. Surface gravity. The surface gravity k of a Killing horizon is defined by the 
formula 

(3.13) d(^giX,X)^ ^ -2kX\ 

where is the one-form metrically dual to X, i.e. X^ = g^^X'^dx^. Two com- 
ments are in order: First, since g{X, X) = on ,yV{X), the differential of g{X, X) 
annihilates T^{X). Now, simple algebra shows that a one-form annihilating a 
null hypersurface is proportional to g{l, •), where ^ is any null vector tangent to 
(those are defined uniquely up to a proportionality factor, see Section 7.1). We 
thus obtain that d{g{X^X)) is proportional to X^] whence (3.13). Next, the name 
"surface gravity" stems from the following: using the Killing equations (3.10) and 

(3.13) one has 

(3.14) X^'W^,X'' = -X^VX^ = kX" . 

Since the left-hand-side of (3.14) is the acceleration of the integral curves of X , the 
equation shows that, in a certain sense, k measures the gravitational field at the 
horizon. 

A key property is that the surface gravity k is constant on bifurcate [240, p. 59] 
Killing horizons. Furthermore, k [209, Theorem 7.1] is constant for all Killing 
horizons, whether bifurcate or not, in space-times satisfying the dominant energy 
condition: this means that 

(3.15) Tp^X^y > for causal future directed vector fields X and Y . 
As an example, consider the Killing vector K of (3.11). We have 

d{g{K, K)) = d{-z^ + t^) = 2(-zrfz + tdt) , 

which equals twice on jV{K)^s. On the other hand, for the Killing vector X of 
(3.12) one obtains 

d{g{X, X)) = 2(< 4- x){dt + dx) , 

which vanishes on each of the Killing horizons {t ~ ~x ,y 7^ 0}. This shows that 
the same null surface can have zero or non-zero values of surface gravity, depending 
upon which Killing vector has been chosen to calculate k. 

The surface gravity of black holes plays an important role in black hole thermo- 
dynamics; see [68] and references therein. 

A Killing horizon .yyo{X) is said to be degenerate, or extreme, if k vanishes 
throughout ,Aq[X); it is called non- degenerate if k has no zeros on jVq{X). Thus, 
the Killing horizons ,jV{K)^s are non-degenerate, while both Killing horizons of X 
given by (3.12) are degenerate. The Schwarzschild black holes have surface gravity 

_ 1 
2m 

So there arc no degenerate black holes within the Schwarzschild family. Theorem 3.5 
below shows that there are no regular, degenerate, static vacuum black holes at all. 

In Kerr space-times we have k = if and only if m = |a|. On the other hand, all 
horizons in the multi-black hole Majumdar-Papapetrou solutions are degenerate. 
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Figure 3.2. A space-like hypersurface M intersecting a Killing 
horizon .yi^{X) in a compact cross-section S. 




Figure 3.3. The general features of the geometry of the orbit- 
space metric on a spacclikc hypersurface intersecting a non- 
degenerate (left) and degenerate (right) Killing horizon, near the 
intersection, visualized by a co-dimension one embedding in Eu- 
clidean space. 

3.4. The orbit-space geometry near Killing horizons. Consider a space-time 
{^,g) with a Killing vector field X. On any set ^ on which X is timelike we can 
introduce coordinates in which X = dt, and the metric may be written as 

(3.16) g = ~V^{dt + e,dx')^ + h,jdx'dx^ , dtV = dtO, ^ dth,j = Q . 

where h = hijdx^dx^ has Ricmannian signature. The metric h is often referred to 
as the orbit-space metric. 

Let M be a spacelike hypersurface in then (3.16) defines a Ricmannian 
metric /i on M n . Assume that X is timelike on a one-sided neighborhood ^ of 
a Killing horizon ,Aq{X), and suppose that M n has a boundary component S 
which forms a compact cross-section of ,A{){X), see Figure 3.2. The vanishing, or 
not, of the surface gravity has a deep impact on the geometry of h near .yfQ{X) [117]: 

(1) Every differentiable such S, included in a degenerate Killing horizon 
.Af){X), corresponds to a complete asymptotic end of {M n '^,h). See 
Figure 3.3.'' 



We are grateful to C. Williams for providing the figure. 
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This remains valid for stationary and axi-symmetric four-dimensional 
configurations without the hypothesis that X is timelike near the hori- 
zon [f39]. 

(2) Every such S included in a smooth Killing horizon ^q{X) on which 

corresponds to a totally geodesic boundary of (M n ^,h), with h being 
smooth up"to~boundary at S. Moreover 

(a) a doubling of (M O ,h) across S leads to a smooth metric on the 
doubled manifold, 

(b) with y^—g{X, X) extending smoothly to — ^y—g{X, X) across S. 

In the Majumdar-Papapetrou solutions of Section 3.9.2, the orbit-space metric h 
as in (3.16) asymptotes to the usual metric on a round cylinder as the event horizon 
is approached. One is therefore tempted to think of degenerate event horizons as 
corresponding to asymptotically cylindrical ends of (A/, h) . 

3.4.1. Average surface gravity. Following [295], near a smooth null hypersurface 
one can introduce Gaussian null coordinates, in which the metric takes the form 

(3.17) g = ripdv'^ + 2dvdr + 2rXadx'^dv + hatdx'^dx'' . 

The hypersurface is given by the equation {r ~ 0}. Let S be any smooth compact 
cross-section of the horizon; then the average surface gravity {k)s is defined as 

(3.18) {k)s = --^^ J^fdfih , 

where d/i/,. is the measure induced by the metric ft, on S*, and \S\ is the volume of 
S. We emphasize that this is defined regardless of whether or not the stationary 
Killing vector is tangent to the null generators of the hypersurface; on the other 
hand, (k)s coincides with k when k is constant and the Killing vector equals d^- 

3.5. Near-horizon geometry. On a degenerate Killing horizon the surface grav- 
ity vanishes, so that the function ip in (3.17) can itself be written as rA, for some 
smooth function A. The vacuum Einstein equations imply (see [295, eq. (2.9)] in 
dimension four and [261, eq. (5.9)] in higher dimensions) 

(3.19) Rab = ^XaXb - Di^aXb) , 

where Rab is the Ricci tensor of hab ■— hab\r=a, and D is the covariant derivative 
thereof, while Xa '■= ha\r=o- The Einstein equations also determine A := A|r=o 
uniquely in terms of ha and hab- 

(3.20) A = h'^" [xaXb - baXb) 

(this equation follows again e.g. from [295, eq. (2.9)] in dimension four, and can 
be checked by a calculation in all higher dimensions). Equations (3.19) have only 
been understood under the supplementary assumptions of staticity [142],^ or axial 
symmetry in space-time dimension four [197,260]: 



^Some partial results with a non-zero cosmological constant have also been proved in [142], 
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Theorem 3.1 ( [142]). Let the space-time dimension be n+l, n>3, suppose that 
a degenerate Killing horizon .vV has a compact cross-section, and that Xa = 9a A 
for some function A (which is necessarily the case in vacuum static space-times). 
Then (3.19) implies Xa = 0, so that hab is Ricci-flat. 

Theorem 3.2 ( [197,260]). In space-time dimension four and in vacuum, suppose 
that a degenerate Killing horizon .jV has a spherical cross- section, and that (./#, g) 
admits a second Killing vector field with periodic orbits. For every connected com- 
ponent jVq of JV there exists a map ip from a neighbourhood of JVq into a degenerate 
Kerr space-time which preserves Xa, hab o,nd A. 

It would be of interest to classify solutions of (3.19) in a useful manner, in 
all dimensions, without any restrictive conditions. Note that the proof that all 
two-dimensional solutions of (3.19) are axi-symmetric would remove the "mean- 
non-degenerate or rotating" restriction from Theorem 3.4 below. 

Theorem 3.2 plays a key role in the proof of uniqueness of stationary axi- 
symmetric degenerate vacuum black holes [139]. Higher-dimensional generalisations 
can be found in [175,211,255,256]. 

In the four-dimensional static case. Theorem 3.1 enforces toroidal topology of 
cross-sections of with a flat hab- On the other hand, in the four-dimensional 
axi-symmetric case. Theorem 3.2 guarantees that the geometry tends to a Kerr one, 
up to second order errors, when the horizon is approached. So, in the degenerate 
case, the vacuum equations impose strong restrictions on the near-horizon geometry. 
This is not the case any more for non-degenerate horizons, at least in the analytic 
setting. 

3.6. Asymptotically flat metrics. In relativity one often needs to consider ini- 
tial data on non-compact manifolds, with natural restrictions on the asymptotic 
geometry. The most commonly studied such examples are asymptotically flat man- 
ifolds, which model isolated gravitational systems. Now, there exist several ways 
of defining asymptotic flatness, all of them roughly equivalent in vacuum. We will 
adapt a Cauchy data point of view, as it appears to be the least restrictive; the 
discussion here will also be relevant for Section 5. 

So, a space-time (^, g) will be said to possess an asymptotically flat end if ./# 
contains a spacelike hypersurface M^xt diffeomorphic to \ B{R), where B{R) is 
a coordinate ball of radius R. An end comes thus equipped with a set of Euclidean 

coordinates {a;*,i = l,...,n}, and one sets r ~ \x\ :~ (X]r=i('^*)^)^^^- ^^'^ then 
assumes that there exists a constant a > such that, in local coordinates on Mext 
obtained from R" \ B{R), the metric h induced by g on M^xt, and the second 
fundamental form K oi M^xt (compare (4.15) below), satisfy the fall-off conditions, 
for some fc > 1, 

(3.21) h,, - ,5,;, = Ofc(r-") , K,j = Ok-i{r-^-") , 

where we write / = Ok{r^) if / satisfies 



(3.22) 



dk,...dkj ^Oir"-') , o<e<k 
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In applications one needs {h, K) to lie in certain weighted Holder or Sobolev space 
defined on Af , with the former better suited for the treatment of the evolution as 
discussed in Section 6.^ 

3.7. Asymptotically flat stationary metrics. For simplicity we assume that 
the space-time is vacuum, though similar results hold in general under appropriate 
conditions on matter fields, see [35, 138] and references therein. 

Along any spacelikc hypersurface M, a Killing vector field X of (^, g) can be 
decomposed as 

X = Nn + Y , 

where Y is tangent to Af , and n is the unit future-directed normal to M . The 
fields N and Y are called "Killing initial data" , or KID for short. The vacuum field 
equations, together with the Killing equations, imply the following set of equations 
on M 

(3.23) DiYj + DjY, = 2NK,j , 

(3.24) R,j{h) + if^Ky - 2i^,feX^, - N-\^yK,j + D.D^N) = , 

where Rij{h) is the Ricci tensor of h. These equations play an important role in 
the gluing constructions described in Section 5.3. 

Under the boundary conditions (3.21), an analysis of these equations provides 
detailed information about the asymptotic behavior of (A^, Y). In particular one can 
prove that if the asymptotic region M^xt is contained in a hypersurface M satisfying 
the requirements of the positive energy theorem (see Section 5.2.1), and if X is 
timelikc along M^xt, then (TV, F*) ~^r~^oo (^°,^*), where the A'^'s are constants 
satisfying {A°)'^ > J2M')'^ [34,138]. Further, in the coordinates of (3.21), 

(3.25) Y' - A' = Okir-") , iV - A" = Ofe(r-") . 

As discussed in more detail in [36], in /i- harmonic coordinates, and in e.g. a maximal 
(i.e., mean curvature zero) time-slicing, the vacuum equations for g form a quasi- 
linear elliptic system with diagonal principal part, with principal symbol identical 
to that of the scalar Laplace operator. It can be shown that, in this "gauge", all 
metric functions have a full asymptotic expansion in terms of powers of In r and 
inverse powers of r. In the new coordinates we can in fact take 

(3.26) a = n-2. 

By inspection of the equations one can further infer that the leading order correc- 
tions in the metric can be written in the Schwarzschild form (3.8). 

3.8. Domains of outer communications, event horizons. A key notion in the 
theory of asymptotically flat black holes is that of the domain of outer communica- 
tions, defined for stationary space-times as follows: For t € K let (/)t[X] : ^ — >■ ^ 
denote the one-parameter group of diffeomorphisms generated by A; we will write 
(j)t for 0t[A] whenever ambiguities are unlikely to occur. Let Mgxt be as in Sec- 
tion 3.6, and assume that X is timelike along M^xt- The exterior region .-#oxt and 



^The analysis of elliptic operators such as the Laplacian on weighted Sobolev spaces was 
initiated by Nirenberg and Walker [303] (see also [23,89,269-272,285-287] as well as [88]). A 
readable treatment of analysis on weighted spaces (not focusing on relativity) can be found in [311]. 
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Figure 3.4. M^xt, ^ext, together with the future and the past of 
^oxt- One has ^ext C /^(^oxt), even though this is not immedi- 
ately apparent from the figure. The domain of outer communica- 
tions is the intersection /+(^oxt)n/~(^oxt), compare Figure 3.5. 



the domain of outer communications ((^cxt)) are then defined as^° 

(3.27) ^cxt := UtMMext) , ((^oxt» = /+(.^oxt) n /"(.^oxt) . 

The black hole region ^ and the black hole event horizon Jif^ are defined as (see 
Figures 3.4 and 3.5) 

(3.28) ^ = ^ \ (.-#ext) , . 

The white hole region W and the white hole event horizon Jif~ are defined as above 
after changing time orientation: 

W = /+(^oxt) , ^"=95^. 

It follows that the boundaries of ((^cxt)) are included in the event horizons. We 
set 

(3.29) §^ = 9((^ext)) n /±(^ext) , s^s+\^s- . 

The sets are achronal boundaries and so. as mentioned in Section 2, they are 
ruled by null geodesies, called generators. 

In general, each asymptotically flat end of ^ determines a different domain of 
outer communications. Although there is considerable freedom in choosing the as- 
ymptotic region Mi^xt giving rise to a particular end, it can be shown that /^(.^cxt), 
and hence ((^cxt)), and §^ . are independent of the choice of M^xt- 

3.9. Uniqueness theorems. It is widely expected that the Kerr metrics provide 
the only stationary, regular, vacuum, four-dimensional black holes. In spite of 
many works on the subject (see [1,81,126,209,224,298,300,336,374] and references 
therein) , the question is far from being settled. See [365] and references therein for 
proposals how to test this idea with astrophysical observations. 

To describe the current state of affairs, some terminology is needed. A Killing 
vector X is said to be complete if its orbits are complete, i.e., for every p S y£ the 
orbit (f)t[X]{j)) of X is defined for all i e M. X is called stationary if it is timelike 
at large distances in the asymptotically flat region. 

A key definition for the uniqueness theory is the following: 



l°See Section 2.3.1 for the definition of /±(f7). 
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Figure 3.5. The hypersurface M from the definition of /+- 
regularity. To avoid ambiguities, we note that ./#oxt is a subset 

of ((^oxt))- 

Definition 3.3. Let {^,g) be a space-time containing an asymptotically flat end 
Mcxt, and let X be a stationary Killing vector field on We will say that 

{^,g,X) is /"^-regular if X is complete, if the domain of outer communications 
((^ext)) is globally hyperbolic, and if ((^ext)) contains a spacelike, connected, 
acausal hypersurface M D M^xt, the closure M of which is a topological manifold 
with boundary, consisting of the union of a compact set and of a finite number of 
asymptotically flat ends, such that the boundary dM := M\M satisfies 

(3.30) dM C S+ 

(see (3.29)), with dM meeting every generator of precisely once; see Figure 3.5. 

Some comments might be helpful. First one requires completeness of the orbits of 
the stationary Killing vector because one needs an action of M on by isometrics. 
Next, one requires global hyperbolicity of the domain of outer communications 
to guarantee its simple connectedness, and to avoid causality violations. Further, 
the existence of a well-behaved spacelike hypersurface gives reasonable control of 
the geometry of ((^oxt))j and is a prerequisite to any elliptic PDEs analysis, as 
is extensively needed for the problem at hand. The existence of compact cross- 
sections of the future event horizon prevents singularities on the future part of 
the boundary of the domain of outer communications, and eventually guarantees 
the smoothness of that boundary. 

The event horizon in a stationary space-time will be said to be rotating if the 
stationary Killing vector is not tangent to the generators of the horizon; it will 
be said mean non- degenerate if < k >gj\77^ (compare (3.18)). The proof of the 
following can be found in [126] in the mean non-degenerate case, and in [139] in 
the degenerate rotating one: 

Theorem 3.4. Let {^,g) he an I'^ -regular, vacuum, analytic, asymptotically flat, 
four- dimensional stationary space-time. If S'~^ is connected and either mean non- 
degenerate or rotating, then ((^oxt)) "is isometric to the domain of outer commu- 
nications of a Kerr space-time. 

Theorem 3.4 finds its roots in work by Carter and Robinson [81,336], with 
further key steps due to Hawking [202] and Sudarsky and Wald [357]. It should be 
emphasized that the hypotheses of connectedness, analyticity, and non-degeneracy 
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ill the case of non-rotating configurations, are higlily unsatisfactory, and one believes 
that they are not needed for the conclusion. 

Recent progress on the connectedness question has been done by Hennig and 
Neugebauer [298] , who excluded two-component configurations under a non-degeneracy 
condition whose meaning remains to be explored; see also [263, 375] for previous 
results. The approach in [298] is based on an observation of Neugebauer and 
Meinel [299] that non-singular two-component configurations, if they exist, must be- 
long to the Neugebauer-Kramer [252] family of solutions. The Neugebauer-Kramer 
metrics are explicitly known, but they contain several free parameters, and the 
metric functions are too complicated to allow a direct analysis of global proper- 
ties. Hennig and Neugebauer exclude most of the solutions by verifying that they 
have negative total ADM mass. For the remaining ones, the authors use a non- 
degeneracy condition introduced by Booth and Fairhurst [54]: a black- hole is said 
to be sub- extremal if any neighborhood of the event horizon contains trapped sur- 
faces. The key of the analysis is a beautiful area inequality of Henning, Ansorg, and 
Cederbaum [207], that on every sub-extremal component of the horizon it holds that 
87r| J| < A, where J is the Komar angular-momentum and A the area of a section. 
Hennig and Neugebauer show that all remaining candidate solutions violate this 
inequality; this is their precise non-existence statement. To tie this with the usual 
theory one would need to show that non-degeneracy in the usual sense implies the 
sub-extremality condition. (The conditions are equivalent for Kerr solutions [54]). 

The analyticity restriction has been removed by Alexakis, lonescu and Klainer- 
man in [1] for ncar-Kerrian configurations, but the general case remains open. The 
key to the approach in [1] is a unique continuation theorem near bifurcate Killing 
horizons proved in [224], which implies the existence of a second Killing vector field, 
say y, in a neighbourhood of the horizon. One then needs to prove that Y extends 
to the whole domain of outer communications. This requires another unique con- 
tinuation theorem [223], which requires specific convexity conditions, and which 
has so far only been shown to apply to ncar-Kerrian configurations. 

Partial results concerning uniqueness of higher dimensional black holes have been 
obtained by Hollands and Yazadjiev [214-216], compare [65,198,199,296]. 

The proof of Theorem 3.4 can be outlined as follows: First, the event horizon in 
a smooth or analytic space-time is a priori only a Lipschitz surface, so the starting 
point of the analysis is provided by a result in [131], that event horizons in regular 
stationary black hole space-times are as differentiable as the differentiability of the 
metric allows. One then shows [202]^^ that either 

a) the stationary Killing vector is tangent to the generators of the event hori- 
zon, or 

b) there exists a second Killing vector defined near the event horizon. 

In case a) one shows that the domain of outer communications contains a maxi- 
mal (mean curvature zero) spacelike hypersurface [145]; to be able to use the result 
from that reference it is useful, first, to extend ((^ext)) using the construction 
in [323]. The construction in [145] assumes a non-degenerate horizon; a proof of 
existence of a maximal hypersurface in the degenerate case is the only element 
needed to remove the "mean-non-degenerate or rotating" restriction of the theo- 
rem. (Alternatively, this restriction could be removed by showing that all solutions 

^-'^Compare [184]; the result, proved by Hawking in space-dimension n = 3 [202,204], has been 
generalised to n > 4 in [212,213,233]. 
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of the near-horizon geometry equations (3.19) are axi-synimetric.) The existence of 
the maximal surface ahows one to estabhsh staticity [357], and one concludes using 
Theorem 3.5 below. 

In case b), analyticity and simple connectedness imply [116] that the isomctry 
group of g) contains a C/(l) factor, with non-empty axis of rotation. A delicate 
argument, which finds its roots in the work of Carter [81], proves that the area 
function 

W ■.= -det{g{Ka,Kb)) , a, 6 =1,2 

where Ka are the stationary and the periodic Killing vector, is strictly positive 
on the domain of outer communications. (This step has been recently simplified 
in [132].) Classical results on group actions on simply connected manifolds [310, 
324] show that the domain of outer communications is diffeomorphic to M x (M'^ \ 
-B(l)), with the action of the isomctry group by translations in the first factor, 
and by rotations around an axis in R'^. The uniformization theorem allows one to 
establish that \/W can be used as the usual polar coordinate p on R^, leading to 
a coordinate system in which the field equations reduce to a harmonic map with 
values in two-dimensional hyperbolic space. The map is singular at the rotation axis 
(compare [121]), with rather delicate singularity structure at points where the event 
horizon meets the axis. A uniqueness theorem for such maps [336,374] achieves the 
proof. 

The analysis above relies heavily on the fact that the domain of outer commu- 
nications is simply connected [146] (compare [188]). 

3.9.1. Static case. Assuming staticity, i.e., stationarity and hypersurface-orthogonality 
of the stationary Killing vector, a more satisfactory result is available in space di- 
mensions less than or equal to seven, and in higher dimensions on manifolds on 
which the Riemannian rigid positive energy theorem holds: non-connected configu- 
rations are excluded, without any a priori restrictions on the gradient of the norm 
of the static Killing vector at event horizons. 

More precisely, we shall say that a manifold M is of positive energy type if there 
are no asymptotically flat complete Riemannian metrics on M with positive scalar 
curvature and vanishing mass except perhaps for a flat one. As made clear in 
Theorem 5.2, this property has been proved so far for all asymptotically flat n- 
dimensional manifolds M of dimension 3 < n < 7 [340] , or under the hypothesis 
that M is spin for any ti > 3, and is expected to be true in general. 

We have the following result, which finds its roots in the work of Israel [236], 
with further simplifications by Robinson [337], and with a significant strengthening 
by Bunting and Masood-ul-Alam [70] ; the proof of the version presented here can 
be found in [117,126], see also [132]: 

Theorem 3.5. Let {.y£',g) be an I'^ -regular, vacuum, static, {n + I) -dimensional 
space-time, n > 3. Let M denote the manifold obtained by doubling the hypersur- 
face M of Definition 3.3 across all non- degenerate components of its boundary and 
smoothly compactifying, in the doubled manifold, all asymptotically flat regions but 
one to a point. If M is of positive energy type, then ((^ext)) is isometric to the 
domain of outer communications of a Schwarzschild space-time. 

Remark 3.6. As a corollary of Theorem 3.5 one obtains non-existence of static, 
regular, vacuum black holes with some components of the horizon degenerate. As 
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observed in [142], if the space-time dimension is four, non-existence follows imme- 
diately from Theorem 3.1 and from simple connectedness of the domain of outer 
communications [146], but this docs not seem to generalize to higher dimensions in 
any obvious way. 

3.9.2. Multi-black hole solutions. In this section wc assume that the space-time 
dimension is four. Space-times containing several black holes seem to be of par- 
ticular interest, but Theorem 3.5 implies, under the conditions spelled out there, 
that no such vacuum solutions exist in the static class. However, the Einstein- 
Maxwell equations admit static solutions with several black holes: the Majumdar- 
Papapetrou (MP) solutions. The metric g and the electromagnetic potential A take 
the form [275,313] 

(3.31) -u-'^dt^ + u'^{dx^ + dy^ + dz'^), A = ±u-^dt, 

with A5U = 0, where A5 is the Laplace operator of the flat metric 5. Standard MP 
black holes are defined by further requiring that 

(3.32) + ^ /fl^, 

for some positive constants /Xj, the electric charges (up to the choice of sign for A 
in (3.31)) carried by the punctures x = Si. Further, the coordinates x'^ of (3.31) 
are required to cover the range M x (R'^ \ {a;}) for a finite set of points di £ M.^, 
i = 1, . . . , I. It has been shown by Hartle and Hawking [200] that standard MP 
space-times can be analytically extended to an electro-vacuum space-time with / 
black hole regions. 

The case / = 1 is the special case m = \q\ of the so-called Reissner- Nordstrom 
metrics, which are the charged, spherically symmetric (connected) generalizations 
of the Schwarzschild black holes with mass m and electric charge q. 

The static /+-regular electro-vacuum black holes are well understood: Indeed, 
the analysis in [81,279,338,355] (compare [118]) leads to: 

Theorem 3.7. Every domain of outer communications in a static, electro-vacuum, 
-regular, analytic black hole space-time without degenerate horizons is isometric 
to a domain of outer communications of a Reissner-Nordstrom black hole. 

The relevance of the standard MP black holes follows now from the following 
result [144]: 

Theorem 3.8. Every domain of outer communications in a static, electro-vacuum, 
-regularblack hole space-time containing degenerate horizons is isometric to a 
domain of outer communications of a standard MP space-time. 

It thus follows that the MP family provides the only static, electro-vacuum, 
analytic, regular black holes with non-connected horizons. 

3.10. Non-zero cosmological constant. A family of black hole solutions with 
non-zero cosmological constant has been discovered by Kottler [251] (compare [50, 
73]). These metrics are also known as the Schwarzschild- anti de Sitter metrics 
(when A < 0) or the Schwarzschild-de Sitter metrics (when A > 0), and take the 
form 

(3.33) g = -e^^^-^dt^ e-2A(.)^^2 ^ ^2^ ^ 
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where k is an Einstein metric on a compact (n— l)-dimensional manifold iV, n > 3. 
Here 

r 

with a = — 2A/n(77,— 1), (3 = i?(fc)/(n— l)(ri, — 2), where R{k) is the scalar curvature 
of the metric k, while m G M is a constant, called the mass of g. The global 
structure of (suitably extended) Kottler space-times has been analyzed in [67, 193] 
in dimension 3 + 1; the results extend to n > 3 dimension. 

4. The Cauchy problem 

The component version of the vacuum Einstein equations with cosmological con- 
stant A (2.8) reads 

(4.1) Ga^+A.g„0 = O, 
where Gap is the Einstein tensor defined as 

(4.2) Gafi := Rap - -^RgaP , 

while Rap is the Ricci tensor and R the scalar curvature. We will refer to those 
equations as the vacuum Einstein equations, regardless of whether or not the cos- 
mological constant vanishes, and in this work we will mostly assume A = 0. Taking 
the trace of (4.1) one obtains 

(«) «=H(ii±i>A. 

n — 1 

where, as elsewhere, n+1 is the dimension of space-time. This leads to the following 
equivalent version of (4.1): 

2A 

(4.4) Ric = -g . 

n — 1 

Thus the Ricci tensor of the metric is proportional to the metric. Pseudo-Riemannian 
manifolds with metrics satisfying Equation (4.4) are called Einstein manifolds in 
the mathematical literature, see e.g. [47]. 

Given a manifold Equation (4.1) or, equivalently. Equation (4.4) forms a 
system of second order partial differential equations for the metric, linear in the 
second derivatives of the metric, with coefficients which are rational functions of 
the ffa/j's, quadratic in the first derivatives of g, again with coefficients rational in g. 
Equations linear in the highest order derivatives are called quasi-linear, hence the 
vacuum Einstein equations constitute a second order system of quasi-linear partial 
differential equations for the metric g. 

In the discussion above we assumed that the manifold ^ has been given. In the 
evolutionary point of view, which we adapt in most of this work, all space-times 
of main interest have topology R x M, where M is an n-dimensional manifold 
carrying initial data. Thus, solutions of the Cauchy problem (as defined precisely 
by Theorem 6.2 below) have topology and differential structure which are deter- 
mined by the initial data. As will be discussed in more detail in Section 6.2, the 
space-times obtained by evolution of the data are sometimes extendible; there is 
then a lot of freedom in the topology of the extended space-time, and we are not 
aware of conditions which would guarantee uniqueness of the extensions. So in the 
evolutionary approach the manifold is best thought of as being given a priori — 
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namely = R x A/, but it should be kept in mind that there is no a priori known 
natural time eoordinate whieh can be constructed by evolutionary methods, and 
which leads to the decomposition ^ = M x M. 

Now, there exist standard classes of partial differential equations which are 
known to have good properties. They are determined by looking at the algebraic 
properties of those terms in the equations which contain derivatives of highest or- 
der, in our case of order two. Inspection of (4.1) shows that this equation does 
not fall in any of the standard classes, such as hyperbolic, parabolic, or elliptic. In 
retrospect this is not surprising, because equations in those classes typically lead 
to unique solutions. On the other hand, given any solution g of the Einstein equa- 
tions (4.4) and any diffeomorphism $, the pull-back metric $*f/ is also a solution of 
(4.4), so whatever uniqueness there might be will hold only up to diffeomorphisms. 
An alternative way of describing this, often found in the physics literature, is the 
following: suppose that we have a matrix g^wix) of functions satisfying (4.1) in 
some coordinate system x^^ . If we perform a coordinate change — >■ y"(x^), then 
the matrix of functions Qapiu) defined as 

Ox^ dx^ 

(4-5) g,,„{x) gaf}{y) = ^^^^(^(y))^^ 

will also solve (4.1), if the x-derivatives there are replaced by y-derivatives. This 
property is known under the name of diffeomorphism invariance, or coordinate 
invariance, of the Einstein equations. Physicists say that "the diffeomorphism 
group is the gauge group of Einstein's theory of gravitation" . 

Somewhat surprisingly, Choquet-Bruhat [178] proved in 1952 that there exists 
a set of hyperbolic equations underlying (4.2). This proceeds by the introduction 
of so-called "harmonic coordinates" , to which we turn our attention in the next 
section. 

4.1. The local evolution problem. 

4.1.1. Wave coordinates. A set of coordinates {y^} is called harmonic if each of 
the functions j/'' satisfies 

(4.6) Dgy^ = , 

where Dg is the d'Alembertian associated with g acting on scalars: 

1 

VI dot .91 

One also refers to these as "wave coordinates". Assuming that (4.6) holds, (4.4) 
can be written as 



(4.7) UJ := tr^Hess / = ^==d^ U\<:\cl g\g^^'' d^f 



n-r 

Here the E^^^'s should be calculated in terms of the g^^'s and their derivatives as in 
(2.3), and the wave operator \3g is as in (4.7). So, in wave coordinates, the Einstein 
equation forms a second-order quasi-linear wave-type system of equations (4.8) for 
the metric functions g°'^ . (This can of course be rewritten as a set of quasi-linear 
equations for the Qaps by algebraic manipulations.) 



28 



PIOTR T. CHRUSCIEL, GREGORY J. GALLOWAY, AND DANIEL POLLACK 



Standard theory of hyperbolic PDEs [173] gives:^^ 
Theorem 4.1. For any initial data 

(4.9) g"^{0,f)eHtt\ dog"^iO,y')eHt,, k > n/2 , 

prescribed on an open subset G C {0} x R" C M x M" there exists a unique solution 
gap gf ^4 defined on an open neighborhood ^ C K x of ^. The set can 
be chosen so that ("^, 9) is globally hyperbolic with Cauchy surface . 

Remark 4.2. The results m [247-249,356] and references therein allow one to reduce 
the differentiability threshold above. 

Equation (4.8) would establish the hyperbolic, evolutionary character of the Ein- 
stein equations, if not for the following problem: Given initial data for an equation 
as in (4.8) there exists a unique solution, at least for some short time. But there is 
a priori no reason to expect that the solution will satisfy (4.6); if it does not, then 
a solution of (4.8) will not solve the Einstein equation. In fact, if we set 

(4.10) := D.y^ , 
then 

1 2 A 

(4.11) i?"^ = -{E"^ - V^A*^ - V^A") + -^^9"^ , 

so that it is precisely the vanishing - or not - of A which decides whether or not a 
solution of (4.8) is a solution of the vacuum Einstein equations. 

This problem has been solved by Choquet-Bruhat [178]. The key observation is 
that (4.11) and the Bianchi identity imply a wave equation for the A"'s. In order 
to see that, recall the twice-contracted Bianchi identity (2.10): 

V„(i?"^'-|5"^)===0. 
Assuming that (4.8) holds, one finds 

= -V„ (V^A'' -f V'A" - V^X-^g"^^ 

= -{n,x^ + R\x") . 

This shows that A" necessarily satisfies the second order hyperbolic system of equa- 
tions 

□gA^ -f i?'^„A" = . 
Now, it is a standard fact in the theory of hyperbolic equations that we will have 

A" = 

on the domain of dependence ^(€?), provided that both A" and its derivatives van- 
ish at i^. To sec how these initial conditions on A" can be ensured, it is convenient 



^■^If k is an integer, then the Sobolev spaces Hj^^^ are defined as spaces of functions which arc in 
L^{K) for any compact set X ,with their distributional derivatives up to order k also in L^(K). In 
the results presented here one can actually allow non-integer fc's, the spaces H^^^ arc then defined 
rather similarly using the Fourier transformation. 
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to assume that y° is the coordinate along the M factor of M x R" , so that the initial 
data surface {0} x is given by the equation = 0. We have 



1 



Vldct.g| ^ 



not n\ \ / 



V^dctgl 

Clearly a necessary condition for the vanishing of Dgy" is that it vanishes at — 0, 
and this allows us to calculate some time derivatives of the metric in terms of space 
ones: 

(4.12) ao(v1drt^.gO") = -a,(v^d^.g-) . 

This implies that the initial data (4.9) for the equation (4.8) cannot be chosen 
arbitrarily if we want both (4.8) and the Einstein equation to be simultaneously 
satisfied. 

Now, there is still freedom left in choosing the wave coordinates. Using this 
freedom, one can show that there is no loss of generality in assuming that on the 
initial hypersurface {y*^ ~ 0} we have 

(4.13) .9°° = -1 , .9"' = , 

and this choice simplifies the algebra considerably. Equation (4.12) determines the 
time derivatives 9o5'^'^|{yO=o} needed in Theorem 4.1, once gij\{yO=o-j and (?o5ij|{yO=o} 
are given. So, from this point of view, the essential initial data for the evolution 
problem arc the space metric 

h ~ g^jdy'dy^ , 

together with its time derivatives. 

It turns out that further constraints arise from the requirement of the vanishing 
of the derivatives of A. Supposing that (4.12) holds at 2/° = — equivalently, 
supposing that A vanishes on {j/*^ = 0}- wc then have 

9,A" = 

on {y^ = 0}. To obtain the vanishing of all derivatives initially it remains to ensure 
that some transverse derivative does. A convenient transverse direction is provided 
by the field n of unit timelike normals to {y° = 0}, and the vanishing of V„A" is 
guaranteed by requiring that 

(4.14) (G^.+Ag^,)n'^ = 0. 

This follows by simple algebra from the equation Eap = and (4.11), 
Gf,i, + Ag^^ = - (^V^A^ + V^A^ - V"Aa5^i. 

using that A^l^o^g = di>^n\yO=o = 0. 

Equations (4.14) are called the Einstein constraint equations, and will be dis- 
cussed in detail in Section 5. 

Summarizing, we have proved: 

Theorem 4.3. Under the hypotheses of Theorem 4-1, suppose that the initial data 
(4.9) satisfy (4.12), (4.13) as well as the constraint equations (4.14). Then the 
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metric given by Theorem 4-1 on a globally hyperbolic set ^ satisfies the vacuum 
Einstein equations. 

4.2. Cauchy data. In Theorem 4.1 we consider initial data given in a single coordi- 
nate patch 1^ C M". This suffices for apphcations such as the Lindblad-Rodnianski 
stabihty theorem discussed in Section 6.5 below, where = R" . But a correct geo- 
metric picture is to start with an rt-dimensional manifold M, and prescribe initial 
data there; the case where M is is thus a special case of this construction. At this 
stage there are two attitudes one may wish to adopt: the first is that M is a subset 
of the space-time ^ — this is essentially what we assumed in Section 4.1. The 
alternative is to consider M as a manifold of its own, equipped with an embedding 

i: M ^ ^ . 

The most convenient approach is to go back and forth between those points of view, 
and this is the strategy that we will follow. 

A vacuum initial data set (M, h, K) is a triple where M is an n-dimensional 
manifold, /i is a Riemannian metric on M, and if is a symmetric two-covariant 
tensor field on M . Further {h, K) are supposed to satisfy the vacuum constraint 
equations that result from (4.14), and which are written explicitly in terms of K 
and h in Section 5.1. Here the tensor field K will eventually become the second 
fundamental form of M in the resulting space-time obtained by evolving the 
initial data. Recall that the second fundamental form of a spacclike hypersurface 
M is defined as 

(4.15) yX G TM K{X, Y) = g{Vxn, Y) , 

where n is the future pointing unit normal to M. The tensor K is often referred to 
as the extrinsic curvature tensor of M in the relativity literature. Specifying K is 
equivalent to prescribing the time-derivatives of the space-part gij of the resulting 
space-time metric g; this can be seen as follows: Suppose, indeed, that a space- 
time {.y^,g) has been constructed (not necessarily vacuum) such that K is the 
extrinsic curvature tensor of M in (^,17). Consider any domain of coordinates 
C M, and construct coordinates in a space-time neighborhood ^ such that 
M n = 1^; those coordinates could be wave coordinates, obtained by solving the 
wave equations (4.6), but this is not necessary at this stage. Since y^ is constant 
on M the one-form dy'^ annihilates TM C T.-#, as does the 1-form g{n, ■). Since 
M has codimension one, it follows that dy*^ must be proportional to g(n, •): 

nady"' = nody" 

on ff. The normalization —1 = g{n,n) = = g^'~'{no)'^ gives 

Uady" = — , dy'^ . 

We then have, by (4.15), 

(4.16) = -^g'^'^(djgai + digaj - dagtj^no . 

This shows that the knowledge of g^i, and dogij at {?/° = 0} allows one to calculate 
Kij. Reciprocally, (4.16) can be rewritten as 

2 

dogij = — Kij + terms determined by the (7/jjy's and their space-derivatives , 
9 ^0 
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SO that the knowledge of the (7p,y's and of the Kij's at = aUows one to calculate 
dogij. Thus, Kij is the geometric counterpart of the dogij's. 

4.3. Solutions global in space. In order to globalize the existence Theorem 4.1 
in space, the key point is to show that two solutions differing only by the values 
9oa\{y'>=o} arc (locally) isometric: so suppose that g and g both solve the vacuum 
Einstein equations in a globally hyperbolic region with the same Cauchy data 
{h, K) on := ^ O M. One can then introduce wave coordinates in a globally 
hyperbolic neighborhood of ^ both for g and g, satisfying (4.13), by solving 

(4.17) Dgy^ = , - , 

with the same initial data for and y^. Transforming both metrics to their 
respective wave-coordinates, one obtains two solutions of the reduced equation (4.8) 
with the same initial data. 

The question then arises whether the resulting metrics will be sufficiently diffcr- 
entiablc to apply the uniqueness part of Theorem 4.1. Now, the metrics obtained 
so far are in a space C^([0, T], i/*), where the Sobolev space iJ^ involves the space- 
derivatives of the metric. The initial data for the solutions or y^^ of (4.17) may 
be chosen to be in H'^'^^ x i/". However, a rough inspection of (4.17) shows that 
the resulting solutions will be only in C^([0, T], i/"), because of the low regular- 
ity of the metric. But then (4.5) implies that the transformed metrics will be in 
C^([0, T], i?"^^), and uniqueness can only be invoked provided thats~\ > n/2 + 1, 
which is one degree of differentiability more than what was required for existence. 
This was the state of affairs for some fifty-five years until the following simple argu- 
ment of Planchon and Rodnianski [320] : To make it clear that the functions y^ are 
considered to be scalars in (4.17), wc shall write y for y^. Commuting derivatives 
with Dg one finds, for metrics satisfying the vacuum Einstein equations, 

□gV„y = V^V^V^y = [V^ V^, Va]y = = . 

Commuting once more one obtains an evolution equation for the field tpa/s '■= 

=0 

where the undcrbraccd term vanishes, for vacuum metrics, by a contracted Bianchi 
identity. So the most offending term in this equation for V'a/Sj involving three 
derivatives of the metric, disappears when the metric is vacuum. Standard theory 
of hyperbohc PDEs shows now that the functions VaV/jy are in C^([0, T], TJ''"^), 
hence y G C^([0, T], i?'*+-^), and the transformed metrics are regular enough to 
invoke uniqueness without having to increase s. 

Suppose, now, that an initial data set {M,h,K) as in Theorem 4.1 is given. 
Covering M by coordinate neighborhoods p € M, one can use Theorem 4.1 to 
construct globally hyperbolic developments {"^p, gp) of (^p, h, K). By the argument 
just given the metrics so obtained will coincide, after performing a suitable coor- 
dinate transformation, wherever simultaneously defined. This allows one to patch 
the (^p,5p)'s together to a globally hyperbolic Lorentzian manifold, with Cauchy 
surface M. Thus: 
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Theorem 4.4. Any vacuum initial data set {M,h,K) of differentiability class 
H"^^ X H" , s > n/2, admits a globally hyperbolic development. 

The solutions are locally unique, in a sense made clear by the proof. The impor- 
tant question of uniqueness in the large will be addressed in Section 6.2. 

4.4. Other hyperbolic reductions. The wave-coordinates approach of Choquct- 
Bruhat, presented above, is the first hyperbolic reduction discovered for the Einstein 
equations. It has been given new life by the Lindblad-Rodnianski stability theorem, 
presented in Section 6.5.2 below. However, one should keep in mind the existence 
of several other such reductions. 

An example is given by the symmetric-hyperbolic first order system of Baum- 
garte, Shapiro, Shibata and Nakamura [30,339,354], known as the BSSN system, 
widely used in numerical general relativity. Another noteworthy example is the 
elliptic- hyperbolic system of [12], in which the elliptic character of some of the 
equations provides increased control of the solution. A notorious problem in numer- 
ical simulations is the lack of constraint preservation, see [217,315] and references 
therein for attempts to improve the situation. The reader is referred to [181, 185] 
for a review of many other possibilities. 

4.5. The characteristic Cauchy problem. Another important systematic con- 
struction of solutions of the vacuum Einstein equations proceeds via a characteristic 
Cauchy problem. In this case the initial data are prescribed on Cauchy hypersur- 
faces which are allowed to be piecewise null. This problem has been considerably 
less studied than the spacelike one described above. We will not go into any details 
here; see [51,71,72,91,92,106,164,326] for further information. 

4.6. Initial-boundary value problems. Numerical simulations necessarily take 
place on a finite grid, which leads to the need of considering initial-boundary value 
problems. In general relativity those are considerably more complicated than the 
Cauchy problem, and much remains to be understood. In pioneering work. Friedrich 
and Nagy [183] constructed a system of equations, equivalent to Einstein's, for a 
set of fields that includes some components of the Weyl tensor, and proved well- 
posedness of an initial-boundary value problem for those equations. It would seem 
that the recent work by Kreiss et al. [253] might lead to a simpler formulation of 
the problem at hand. 

5. Initial data sets 

We now turn our attention to an analysis of the constraint equations, returning 
to the evolution problem in Section 6. 

An essential part of the mathematical analysis of the Einstein field equations 
of general relativity is the rigorous formulation of the Cauchy problem, which is a 
means to describe solutions of a dynamical theory via the specification of initial data 
and the evolution of that data. In this section we will be mainly concerned with 
the initial data sets for the Cauchy problem. As already explained in Section 4.1.1, 
those initial data sets have to satisfy the relativistic constraint equations (4.14). 
This leads to the following questions: What are the sets of allowable initial data? 
Is it possible to parameterize them in a useful way? What global properties of the 
space-time can be seen in the initial data sets? How does one engineer initial data 
so that the associated space-time has some specific properties? 
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5.1. The constraint equations. As explained in Section 4.2, an initial data set 
for a vacuum space-time consists of an n-dimensional manifold M together with a 
Riemannian metric h and a symmetric tensor K. In the non-vacuum case we also 
have a collection of non-gravitational fields which we collectively label F (usually 
these are sections of a bundle over M). We have already seen the relativistic 
vacuum constraint equations expressed as the vanishing of the normal components 
of the Einstein equations (4.14). Now, if h is the metric induced on a spacclike 
hypersurface in a Lorentzian manifold, it has its own curvature tensor R^jki- If we 
denote by Kij the second fundamental form of M in and by ^^jki the space-time 
curvature tensor, the Gauss-Codazzi equations provide the following relationships: 



Here n is the timelike normal to the hypersurface, and we are using a coordinate 
system in which the 9i's are tangent to the hypersurface M. 

Contractions of (5.1)-(5.2) and simple algebra allow one to reexpress (4.14) in 
the following form, where we have now allowed for the additional presence of non- 
gravitational fields: 



where R{h) is the scalar curvature of the metric h, J is the momentum density 
of the non-gravitational fields, p is the energy density,^'^ and C{J-,h) denotes the 
set of additional constraints that might come from the non-gravitational part of 
the theory. The first of these equations is known as the momentum constraint and 
is a vector field equation on M. The second, a scalar equation, is referred to as 
the scalar, or Hamiltonian, constraint, while the last are collectively labeled the 
non-gravitational constraints. These are what we shall henceforth call the Einstein 
constraint equations, or simply the constraint equations if ambiguities are unlikely 
to occur. 

As an example, for the Einstein-Maxwell theory in 3+1 dimensions, the non- 
gravitational fields consist of the electric and magnetic vector fields E and B. In 
this case we have p = ^{\E\1 + \B\fJ, J = [E x B)h, and we have the extra 
(non-gravitational) constraints div^ E ~ and divh B — 0. 

Equations (5.3)-(5.5) form an underdetermined system of partial differential 
equations. In the classical vacuum setting of n = 3 dimensions, these are locally 
four equations for the twelve unknowns given by the components of the symmetric 
tensors h and K. This section will focus primarily on the vacuum case with a zero 
cosmological constant. However, we will allow arbitrary values of A in Section 5.3.5. 

The most successful approach so far for studying the existence and uniqueness 
of solutions to (5.3)-(5.5) is through the conformal method of Lichnerowicz [264], 
Choquet-Bruhat and York [100]. The idea is to introduce a set of unconstrained 
"conformal data" , which are freely chosen, and find {h, K) by solving a system 



If T is the stress-energy tensor of the non-gravitational fields, and n denotes the unit timchke 
normal to a hypersurface M embedded in a space-time, with induced data {M,h,K,T), then 
J = —T{n,-) and p = T{n,n). However, in terms of the initial data sot itself we shall regard 
(5.3)-(5.4) as the definitions of the quantities J and p. 



(5.1) 
(5.2) 



R'jki = ^'jki + K\Kjk - K\Kji 



(5.3) 
(5.4) 
(5.5) 



dWK - d{trK) 
Rih) - 2A - \K\l + {trKf 
C{T, h) 



Sir J , 
IGnp , 
0, 
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of determined partial differential equations. In the vacuum case with vanishing 
cosmological constant [100], the free conformal data consist of a manifold M, a 
Riemannian metric h on M, a trace-free symmetric tensor a, and the mean curva- 
ture function t. The initial data {h, K) defined as 

(5.6) h = cp'ih, q=^, 

n — 2 

(5.7) K = (j)-^{d + VW) + -(t)'ih, 

n 

where cj) is positive, will then solve (5.3)-(5.4) if and only if the function (j) and the 
vector field W solve the equations 

~ 77 — 1 

(5.8) divr (VW + ct) = (f)i+^DT , 

n 

(5.9) - —±—R{h)^ + + VW\lr'-^ - — r20«+i = . 

q[n — 1) q[n — 1) " qn 

We use the symbol D to denote the covariant derivative of /i; T) is the conformal 
Killing operator. 

(5.10) VWab = DaWb + bbWa - -habDcW^ ■ 

n 

Vector fields W annihilated by V are called conformal Killing vector fields, and are 
characterized by the fact that they generate (perhaps local) conformal diffeomor- 
phisms of (M, h). The semi-linear scalar equation (5.9) is often referred to as the 
Lichnerowicz equation. 

Equations (5.8)-(5.9) form a determined system of equations for the (n + 1) 
functions {(f>,W). The operator div,-(2? •) is a linear, formally self-adjoint, elliptic 
operator on vector fields. What makes the study of the system (5.8)-(5.9) difficult 
in general is the nonlinear coupling between the two equations. 

The explicit choice of (5.6)-(5.7) is motivated by the two identities (for h = (j)'^h) 

(5.11) R{h) - (g(n - 1) A,,0 - Rih)<f>) , 

where q = ;;732' "^hich is the unique exponent that does not lead to supplementary 
terms in (5.11), and 

(5.12) Dlir'^Bab) = ^-''-^D-^Bab 

which holds for any trace-free tensor B. Equation (5.11) is the well known identity 
relating the scalar curvatures of two conformally related metrics. 

In the space-time evolution (^, g) of the initial data set (M, h, K), the function 
T = tr^K is the mean curvature of the hypersurface M C The assumption that 
the mean curvature function r is constant on M significantly simplifies the analysis 
of the vacuum constraint equations because it decouples equations (5.8) and (5.9). 
One can then attempt to solve (5.8) for W, and then solve the Lichnerowicz equation 
(5.9).^ 

Existence and uniqueness of solutions of this problem for constant mean curva- 
ture ("CMC") data has been studied extensively. For compact manifolds this was 
exhaustively analysed by Isenberg [225], building upon a large amount of previ- 
ous work [100,264,306,380]; the proof was simplified by Maxwell in [280]. If we let 
3^([ft-]) denote the Yamabe invariant of the conformal class [h] of metrics determined 
by h (see [257]), the result reads as follows: 
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Theorem 5.1 ( [225]). Consider a smooth conformal initial data set {h,a,T) on 
a compact manifold M , with constant t . Then there always exists a solution W of 
(5.8). Setting a = VW + a, the existence, or not, of a positive solution (jj of the 
Lichnerowicz equation is shown in Table 5.1. 





cr = 0,T = 


(J = 0,T^O 


cr ^ 0,r = 


cr^O,T^O 


y{[h]) < 


No 


Yes 


No 


Yes 


ym) = o 


Yes 


No 


No 


Yes 


y{[h]) > 


No 


No 


Yes 


Yes 



Table 5.1. Existence of solutions in the conformal method for 
CMC data on compact manifolds. 



More recently, work has been done on analyzing these equations for metrics of low 
differentiability [86,280]; this was motivated in part by recent work on the evolution 
problem for "rough initial data" [247-249,356]. Exterior boundary value problems 
for the constraint equations, with nonlinear boundary conditions motivated by black 
holes, were considered in [157,281]. 

The conformal method easily extends to CMC constraint equations for some 
non-vacuum initial data, e.g. the Einstein-Maxwell system [225] where one obtains 
results very similar to those of Theorem 5.1. However, other important examples, 
such as the Einstein-scalar field system [95-97, 205] , require more effort and are not 
as fully understood. 

Conformal data close to being CMC (e.g. via a smallness assumption on [Vr]) 
are usually referred to as "near-CMC" . Classes of near-CMC conformal data solu- 
tions have been constructed [2,94,230,231] and there is at least one example of a 
non-existence theorem [234] for a class of near-CMC conformal data. However, due 
to the non-linear coupling in the system (5.8)-(5.9), the question of existence for un- 
restricted choices of the mean curvature r appears to be significantly more difficult, 
and until recently all results assumed strong restrictions on the gradient of r. The 
first general result in this context is due to Hoist, Nagy, and Tsogtgerel [218,219], 
who construct solutions with freely specified mean curvature in the presence of mat- 
ter. In [284], Maxwell provides a sufficient condition, with no restrictions on the 
mean curvature, for the conformal method to generate solutions to the vacuum con- 
straint equations on compact manifolds. As an application. Maxwell demonstrates 
the existence of a large class of solutions to the vacuum constraint equations with 
freely specified mean curvature. These results together represent a significant ad- 
vance in our understanding of how the conformal method may be used to generate 
solutions of the vacuum constraint equations. However the existence question for 
generic classes of large conformal data remains wide open; compare [283] for some 
new results. 

The analysis of the conformal constraint equations (5.8)-(5.9) discussed above 
proceeds either via the method of sub- and super-solutions (which is a barrier 
argument exploiting the maximum principle), or a perturbation or fixed point 
method. In [205] Hebey, Pacard and Pollack used the mountain pass lemma to 
analyse Lichnerowicz-type equations arising in certain cases of the Einstein-scalar 
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field system. Such arguments may conceivably prove useful in studying (5.8)-(5.9) 
for general r's. 

A natural question is whether the set of solutions to the constraint equations 
forms a manifold. This was first considered by Fisher and Marsden [177], who pro- 
vided a Frcchet manifold structure; Banach manifold structures have been obtained 
in [129], and a Hilbert manifold structure (for asymptotically flat initial data sets) 
in [27]. 

In [29] the reader will find a presentation of alternative approaches to construct- 
ing solutions of the constraints, covering work done up to 2003. 

5.1.1. The constraint equations on asymptotically flat manifolds. Recall that asymp- 
totically flat initial data sets have been defined in Section 3.6. There is a large 
number of well-established results concerning the existence of CMC and "near 
CMC" solutions of the Einstein constraint equations on asymptotically fiat man- 
ifolds [78,88,89,98,100,281,282]. We mention here only that there is a precise 
conformally invariant criterion sufficient to prove the existence of a positive so- 
lution to the Lichnerowicz equation relative to a given set of asymptotically fiat 
conformal data [78,281]". 

5.2. Mass inequalities. Among the deepest results in mathematical general rel- 
ativity are the global mass inequalities for asymptotically fiat manifolds. Those 
have been discussed extensively in the existing literature [57,58,63,257,340,342], 
and therefore will only be given the minimum amount of attention, as needed for 
the remaining purposes of this work. 

5.2.1. The Positive Mass Theorem. Using the coordinate system of (3.21), one 
defines the Arnowitt-Deser-Misner [15] mass of (M, h) of an asymptotically flat 
end as 

Here Sr is the coordinate sphere at radius r and dSi = di\d^, and d^ is the 
Riemannian volume form of h. The factor IGtt is a matter of convention and is 
natural in space-dimension three. The integral converges to a finite, coordinate- 
independent limit if, for some a > -^^y^, 

(5.14) \h,j - 5,j\ < cr-'', \dh\ < cr-"-^ and Rih) G L\M,,t), 

with those conditions being essentially optimal [23,112]. 

For time-symmetric initial data the vacuum constraint equations (5.3)-(5.4) re- 
duce to the condition that the metric h be scalar-fiat, i.e. R{h) = 0. On the other 
hand, if one considers time-symmetric data for a non- vacuum space-time, then from 
(5.4) we see that the scalar curvature is twice the energy density of the matter fields. 
The non-negativity of R, assuming a vanishing cosmological constant A, is then a 
consequence of the dominant energy condition for initial data (which follows from 

(3.15) ), 

(5.15) P>\J\h, 



^^Maxwell [281] uncovered an error in Cantor's definition of the invariant [78] and provided 
the correct definition. 
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where p and J are defined in (5.3)-(5.4); see also Section 7.5. One checks that the 
dominant energy condition (3.15) holds on ^ if and only if (5.15) holds relative 
to each spacelikc hypcrsurfacc in ^ . 

Now, the ADM mass is thought to represent the total mass of the system as 
viewed on Af , which contains contributions from the matter fields, the gravitational 
field, as well as their binding energy. The long-standing question of its positivity 
was resolved by Schocn and Yau [344] in dimension three, and is now known as the 
Positive Mass Theorem: 

Theorem 5.2. Let {M,h) be an asymptotically flat Riemannian manifold with 
nonnegative scalar curvature. Suppose that either M is spin, or the dimension 
n <1 , or that h is conformally flat. Then the total ADM mass m satisfies m > 0, 
with equality if and only if{M,h) is isometric to Euclidean space (]R",(5). 

As remarked in the introduction, this theorem, and its generalizations, stands 
as one of the cornerstones of mathematical relativity. Accessible introductions to 
the positive mass theorem may be found in [58, 119,257,340]. The restriction on 
the dimension arises from the use of area minimizing hypersurfaces [340] , which are 
known to sometimes possess singularities in higher dimensions, and it is expected 
that positivity is true in all dimensions. The Positive Mass Theorem was proven in 
all dimensions for conformally flat manifolds by Schoen and Yau [349] by a different 
argument, and in all dimensions for spin manifolds by Witten [378] (see also [23, 
85,314]). The result was generalized in [346,347] (compare [342]) to asymptotically 
flat initial data sets {M,h,K,T) satisfying the dominant energy condition (5.15). 

5.2.2. Riemannian Penrose Inequality. An important generalization of the Positive 
Mass Theorem is given by the Riemannian Penrose Inequality. 

Theorem 5.3. Let (M, h) he a complete, smooth, asymptotically flat 3-manifold 
with nonnegative scalar curvature with total mass m and which has an outermost 
minimal surface Sq of area Aq . Then 

(5.16) ^>/A, 

with equality if and only if (M, h) is isometric to the Schwarzschild metric (M'^ \ 
{0}, (1 + 2pj)^'^) outside their respective outermost minimal surfaces. 

Theorem 5.3 was first proved by Huisken and Ilmancn [222] under the restric- 
tion that Aq is connected, or assuming instead that Aq is the area of the largest 
connected component of Eg. The version above, with a proof that uses completely 
different methods, is due to Bray [55]. The proofs are beautiful applications of 
geometric flows to a fundamental problem in relativity. A number of accessible re- 
views has been written on these important results, to which we refer the interested 
reader [56-58,63,276]. A generalization of Theorem 5.3 to dimensions n < 7 has 
been established in [61]. 

One expects that some form of (5.16) holds for general relativistic initial data sets 
(h, K) satisfying the dominant energy condition. A suggestion how one could prove 
this has been put forward by Bray and Khuri in [59,60], compare [64,79,179,277]. 
An inequality in the spirit of (5.16), but involving some further geometric constants, 
has been proved by Herzlich [208] in the Riemannian case. 
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5.2.3. Quasi-local mass. In the context of asymptotically flat space-times, there are 
well-defined global notions of mass and energy, and these are central to the cele- 
brated positive mass theorem discussed in Section 5.2.1. One would, however, like 
to have a well-defined useful local notion of mass or energy, with natural properties 
- e.g., monotonicity - that one has in other physical theories. Such a definition 
has been elusive despite a great deal of effort by many people and this remains 
an important open problem. We refer the reader to the Living Reviews article 
by Szabados [359] for a survey, and note that there have been interesting recent 
mathematical developments in the area [62,288,353,371-373], not described in the 
currently available version of [359] . 

5.3. Applications of gluing techniques. Over the past 25 years, "gluing tech- 
niques" have become a standard tool in geometric analysis. Since the construc- 
tion by Taubes of self-dual Yang-Mills connections on four- manifolds [364], which 
played a crucial role in Donaldson's construction of exotic smooth structures in 
four-dimensions [162], gluing has been applied in important ways across a very 
wide range of areas. What gluing typically refers to is a construction in which 
solutions of a nonlinear partial differential equation or system, which correspond to 
some geometric quantity of interest, e.g. self-dual connections, are fused together to 
create new solutions. This is done by a mix of geometry and analysis in which one 
ultimately studies the linearization of the relevant PDEs, and in most cases one has 
to overcome analytic degeneracy introduced in the gluing procedure. Thus, from an 
analytic point of view, gluing should be regarded as a singular perturbation method. 
Part of the usefulness of the technique lies in the fact that, away from the small 
set about which one fuses the two solutions, the new solution is very close to the 
original ones. The fact that the original solutions are usually not exactly preserved 
is a reflection of the fact that the relevant equations satisfy a unique continuation 
property: any two solutions which agree on an open set must agree everywhere. 
This is a well known property for, say, a scalar semi-linear elliptic equation. The 
constraint equations, however, are underdetermined. As explained below, Corvino 
and Schoen have introduced gluing techniques which exploit this and show the fail- 
ure (to a high degree) of unique continuation. This leads to localized gluings of 
initial data sets, which has proven very valuable. 

5.3.1. The linearized constraint equations and KIDs. The starting point of gluing 
constructions for the constraint equations is the linearization of these equations 
about a given solution {M,h,K). Wc let V*/-^ j^-^ denote the adjoint of the lin- 
earization of the constraint equations at this solution. Viewed as an operator acting 
on a scalar function N and a vector field Y, VT/^ takes the explicit form [128] 



Now this operator docs not. on first inspection, appear to be very "user friendly". 
However, our immediate concern is solely with its kernel, and the pairs (TV, Y) which 
lie in its kernel have a very straightforward geometric and physical characterization. 
In particular, let 11 be an open subset of M. By definition, the set of "KIDs" on 



/ 2(V(,y,) - W'Yig,, - K,jN + tvK Ng,,) 



\ 



(5.17) Vl^^K^iN,Y) 



W'YiK,, - 2K\,^^)Yi + Kh^qY'g^3 
-ANg,, + V,V,iV + (VP^ipffy - ViK,,)Y' 
\ -iVRic {g),j + 2NK\Kji - 27V(tr K)Kij 
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fl, denoted J^{^), is the set of all solutions of the equation 
(5.18) 7'(V^)|JiV,r) = 0. 

Such a solution {N,Y), if nontrivial, generates a space-time Killing vector field in 
the domain of dependence of {il, K\q) [291]; compare Section 3.6. 

From a geometric point of view one expects that solutions with symmetries 
should be rare. This was made rigorous in [37], where it is shown that the generic 
behaviour among solutions of the constraint equations is the absence of KIDs on 
any open set. On the other hand, one should note that essentially every explicit 
solution has symmetries. In particular, both the flat initial data for Minkowski 
space, and the initial data representing the constant time slices of Schwarzschild 
have KIDs. 

5.3.2. Corvino's result. As we have already pointed out, the Einstein constraint 
equations form an underdetermined system of equations, and as such, it is unrea- 
sonable to expect that they (or their linearizations) should satisfy the unique contin- 
uation property. In 2000, Corvino established a gluing result for asymptotically flat 
metrics with zero scalar curvature which dramatically illustrated this point [148]. 
In the special case when one considers initial data with vanishing second funda- 
mental form K = 0. the momentum constraint equation (5.3) becomes trivial and 
the Hamiltonian constraint equation (5.3) reduces to simply R(h) = 0, i.e. a scalar 
flat metric. Such initial data sets are referred to as "time-symmetric" because the 
space-time obtained by evolving them possesses a time-reversing isometry which 
leaves the initial data surface flxed. Beyond Euclidean space itself, the constant 
time slices of the Schwarzschild space-time form the most basic examples of asymp- 
totically flat, scalar flat manifolds. One long-standing open problem [26,350] in the 
fleld had been whether there exist scalar flat metrics on M" which are not globally 
spherically symmetric but which are spherically symmetric in a neighborhood of 
infinity and hence, by Birkhoff's theorem, Schwarzschild there. 

Corvino resolved this by showing that he could deform any asymptotically fiat, 
scalar flat metric to one which is exactly Schwarzschild outside of a compact set. 

Theorem 5.4 ( [148]). Let (M,h) be a smooth Riemannian manifold with zero 
scalar curvature containing an asymptotically flat end M^xt = {l^;] > r > 0}. Then 
there is a R > r and a smooth metric h on M with zero scalar curvature such that h 
is equal to h in M\Mcxt o,nd h coincides on {\x\ > R} with the metric induced on 
a standard time- symmetric slice in the Schwarzschild solution. Moreover the mass 
of h can be made arbitrarily close to that of h by choosing R sufficiently large. 

Underlying this result is a gluing construction where the deformation has com- 
pact support. The ability to do this is a reflection of the underdetermined nature 
of the constraint equations. In this setting, since K = 0, the operator takes a much 
simpler form, as a two-covariant tensor valued operator acting on a scalar function 
u by 

V*u = —{Ahu)h + Hess/iU — uRic{h) . 

An elementary illustration of how an underdetermined system can lead to compactly 
supported solutions is given by the construction of compactly supported transverse- 
traceless tensors on R'^ in Appendix B of [149] (see also [33, 158]). 

An additional challenge in proving Theorem 5.4 is the presence of KIDs on the 
standard slice of the Schwarzschild solution. If the original metric had ADM mass 
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m{h), a naive guess could be that the best fittmg Schwarzschild solution would 
be the one with precisely the same mass. However the mass, and the coordinates 
of the center of mass, are in onc-to-onc correspondence with obstructions arising 
from KIDs. To compensate for this co-kernel in the linearized problem, Corvino 
uses these (n + 1 in dimension n) degrees of freedom as efFcctive parameters in the 
geometric construction. The final solution can be chosen to have its ADM mass 
arbitrarily close to the initial one. 

The method uncovered in Corvino's thesis has been applied and extended in a 
number of important ways. The "asymptotic simplicity" model for isolated gravi- 
tational systems proposed by Penrose [316] has been very influential. This model 
assumes existence of smooth conformal completions to study global properties of 
asymptotically flat space-times. The question of existence of such vacuum space- 
times was open until Chrusciel and Delay [127], and subsequently Corvino [149], 
used this type of gluing construction to demonstrate the existence of infinite di- 
mensional families of vacuum initial data sets which evolve to asymptotically sim- 
ple space-times. The extension of the gluing method to non-timc-symmctric data 
was done in [128,150]. This allowed for the construction of space-times which are 
exactly Kerr outside of a compact set, as well as showing that one can specify other 
types of useful asymptotic behavior. 

5.3.3. Conformal gluing. In [227], Isenberg, Mazzeo and Pollack developed a gluing 
construction for initial data sets satisfying certain natural non-degeneracy assump- 
tions. The perspective taken there was to work within the conformal method, and 
thereby establish a gluing theorem for solutions of the determined system of PDEs 
given by (5.8) and (5.9). This was initially done only within the setting of constant 
mean curvature initial data sets and in dimension n = 3 (the method was extended 
to all higher dimensions in [226]). The construction of [227] allowed one to combine 
initial data sets by taking a connected sum of their underlying manifolds, to add 
wormholes (by performing codimension 3 surgery on the underlying, connected, 3- 
manifold) to a given initial data set, and to replace arbitrary small neighborhoods 
of points in an initial data set with asymptotically hyperbolic ends. 

In [228] this gluing construction was extended to only require that the mean 
curvature be constant in a small neighborhood of the point about which one wanted 
to perform a connected sum. This extension enabled the authors to show that one 
can replace an arbitrary small neighborhood of a generic point in any initial data 
set with an asymptotically flat end. Since it is easy to see that CMC solutions of 
the vacuum constraint equations exist on any compact manifold [377], this leads 
to the following result which asserts that there arc no topological obstructions to 
asymptotically fiat solutions of the constraint equations. 

Theorem 5.5 ( [228]). Let M be any closed n-dimensional manifold, and p G M. 
Then M \ {p} admits an asymptotically flat initial data set satisfying the vacuum 
constraint equations. 

5.3.4. Initial data engineering. The gluing constructions of [227] and [228] are per- 
formed using a determined elliptic system provided by the conformal method, which 
necessarily leads to a global deformation of the initial data set, small away from the 
gluing site. Now, the ability of the Corvino gluing technique to establish compactly 
supported deformations invited the question of whether these conformal gluings 
could be localized. This was answered in the affirmative in [128] for CMC initial 
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data under the additional, generically satisfied [37], assumption that there are no 
KIDs in a neighborhood of the gluing site. 

In [135,136], this was substantially improved upon by combining the gluing 
construction of [227] together with the Corvino gluing technique of [127,148], to 
obtain a localized gluing construction in which the only assumption is the absence of 
KIDs near points. For a given n-nianifold M (which may or may not be connected) 
and two points Pa G M, a = 1,2, we let M denote the manifold obtained by 
replacing small geodesic balls around these points by a neck S*""^ x /. When M is 
connected this corresponds to performing codimension n surgery on the manifold. 
When the points Pa lie in different connected components of M, this corresponds 
to taking the connected sum of those components. 

Theorem 5.6 ( [135,136]). Let (il/, h,K) be a smooth vacuum initial data set, with 
M not necessarily connected, and consider two open sets fla C M, a ~ 1,2, with 
compact closure and smooth boundary such that 

the set of KIDs, J(f{^la), is trivial. 

Then for all pa € fia. e > 0, and fc g N there exists a smooth vacuum initial data 
set {M , h{e) , K [e)) on the glued manifold M such that {h{e),K{e)) is e-close to 
{h, K) in a x topology away from B{pi, e) U B{p2, e). Moreover (/i(e), K{e)) 
coincides with {h, K) away from ili U 172- 

This result is sharp in the following sense: first note that, by the positive mass 
theorem, initial data for Minkowski space-time cannot locally be glued to anything 
which is non-singular and vacuum. This meshes with the fact that for Minkowskian 
initial data, we have J^(r2) ^ {0} for any open set J7. Next, recall that by the results 
in [37], the no-KIDs hypothesis in Theorem 5.6 is generically satisfied. Thus, the 
result can be interpreted as the statement that for generic vacuum initial data sets 
the local gluing can be performed around arbitrarily chosen points Pa- In particular 
the collection of initial data with generic regions fia satisfying the hypotheses of 
Theorem 5.6 is not empty. 

The proof of Theorem 5.6 is a mixture of gluing techniques developed in [226,227] 
and those of [128, 148, 150]. In fact, the proof proceeds initially via a generalization 
of the analysis in [227] to compact manifolds with boundary. In order to have CMC 
initial data near the gluing points, which the analysis based on [227] requires, one 
makes use of the work of Bartnik [24] on the plateau problem for prescribed mean 
curvature spacelike hypersurfaces in a Lorentzian manifold. 

Arguments in the spirit of those of the proof of Theorem 5.6 lead to the construc- 
tion of many-body initial data [124,125]: starting from initial data for N gravitating 
isolated systems, one can construct a new initial data set which comprises isomet- 
rically compact subsets of each of the original systems, as large as desired, in a 
distant configuration. 

An application of the gluing techniques concerns the question of the existence of 
CMC slices in space-times with compact Cauchy surfaces. In [25], Bartnik showed 
that there exist maximally extended, globally hyperbolic solutions of the Einstein 
equations with dust which admit no CMC slices. Later, Eardley and Witt (un- 
published) proposed a scheme for showing that similar vacuum solutions exist, but 
their argument was incomplete. It turns out that these ideas can be implemented 
using Theorem 5.6, which leads to: 
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Corollary 5.7. [135,136] There exist maximal globally hyperbolic vacuum space- 
times with compact Cauchy surfaces which contain no compact spacelike hypersur- 
faces with constant mean curvature. 

Compact Cauchy surfaces with constant mean curvature are useful objects, as 
the existence of one such surface gives rise to a unique fohation by such surfaces [66] , 
and hence a canonical choice of time function (often referred to as CMC or York 
time). Foliations by CMC Cauchy surfaces have also been extensively used in 
numerical analysis to explore the nature of cosmological singularities. Thus the 
demonstration that there exist space-times with no such surfaces has a negative 
impact on such studies. 

One natural question is the extent to which space-times with no CMC slices 
are common among solutions to the vacuum Einstein equations with a fixed spatial 
topology. It is expected that the examples constructed in [135,136] are not isolated. 
In general, there is a great deal of flexibility (in the way of free parameters) in the 
local gluing construction. This can be used to produce one parameter families of 
distinct sets of vacuum initial data which lead to space-times as in Corollary 5.7. 
What is less obvious is how to prove that all members of these families give rise to 
distinct maximally extended, globally hyperbolic vacuum space-times. 

A deeper question is whether a sequence of space-times which admit constant 
mean curvature Cauchy surfaces may converge, in a strong topology, to one which 
admits no such Cauchy surface. (See [22,25, 191] for general criteria leading to the 
existence of CMC Cauchy surfaces.) 

5.3.5. Non-zero cosmological constant. Gluing constructions have also been carried 
out with a non-zero cosmological constant [130, 140, 141]. One aim is to construct 
space-times which coincide, in the asymptotic region, with the corresponding black 
hole models. In such space-times one has complete control of the geometry in the 
domain of dependence of the asymptotic region, described there by the Kottler met- 
rics (3.33). For time-symmetric slices of these space-times, the constraint equations 
reduce to the equation for constant scalar curvature R = 2A. Gluing constructions 
have been previously carried out in this context, especially in the case of A > 0, 
but in [130, 140, 141] the emphasis is on gluing with compact support, in the spirit 
of Corvino's thesis and its extensions already discussed. 

The time-symmetric slices of the A > Kottler space-times provide "Delaunay" 
metrics (see [141] and references therein), and the main result of [140,141] is the 
construction of large families of metrics with exactly Delaunay ends. When A < 
the focus is on asymptotically hyperbolic metrics with constant negative scalar 
curvature. With hindsight, within the family of Kottler metrics with A £ M (with 
A = corresponding to the Schwarzschild metric), the gluing in the A > setting 
is technically easiest, while that with A < is the most difficult. This is due to the 
fact that for A > one deals with one linearized operator with a one-dimensional 
kernel; in the case A = the kernel is (n-l-l)-dimensional; while for A < one needs 
to consider a one-parameter family of operators with (n + l)-dimensional kernels. 

6. Evolution 

In Section 4 we have seen that solutions of the vacuum Einstein equations can 
be constructed by solving a Cauchy problem. It is then of interest to inquire 
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about the global properties of the resulting space-times. As an illustration of 
what can happen, we start with a short description of a solution of the vacuum 
Einstein equations with interesting dynamical properties, the Taub - Newman Unti 
Tamburino (Taub-NUT) metric. 

6.1. Taub— NUT space-times. The Taub-NUT metrics provide interesting ex- 
amples of pathological behavior: incomplete geodesies within compact sets, closed 
causal curves, inequivalent extensions of maximal globally hyperbolic regions (to be 
defined shortly) and inequivalent conformal boundary completions at infinity. For 
all these reasons they are a must in any introductory discussion of general relativity. 

The Taub-NUT metrics [301,363] are solutions of the vacuum Einstein equations 
on space-time manifolds of the form 

.^i := / X , 

where / is an interval. They take the form [290] 

(6.1) ^U-^de + {2lfUa\ + (<2 + e){al + af) , 

(6-2) V{t) = -\ + ^^^^. 

Here £ and m are real numbers with £ > 0. Further, the one-forms (Xi, ai and (73 
form a basis for the set of left-invariant one- forms on SU{2) S^: If 

is3 -.S^^W^ 

is the standard embedding of S"^ into W^, then one can take 

(6.3) (Ti = 2ig3{x dw — w dx + y dz ~ z dy) , 

(72 = 2ig3{z dx — X dz + y dw — w dy) , 
(73 = 2ig3{x dy — y dx + z dw — w dz) . 

The metric (6.1) is invariant under a left SU (2)-action, and a further U (1) action 
consisting of right-rotations of (72 and (73 amongst themselves, so that the connected 
component of the identity of the group of isometries of g is SU{2) x U{1). 

The function U always has two zeros, 

_ (t+-t)(t-t_) 

where 



t± := TO ± \/m^ + £^ . 

It follows that / has to be chosen so that t± ^ I. The space-time {^(t-.t+)i 9) will 
be referred to as the Taub space-time [363]. It is not very difficult to show that the 
Taub space-times exhaust the collection of maximal globally hyperbolic vacuum 
space-times evolving from SU{2) x J7(l)"invariant Cauchy data on . 

The significance of these metrics to strong cosmic censorship, which we are about 
to discuss, stems from the following observation (see [120,133,290]): 

Theorem 6.1. Consider the Taub space-time {■^(t-,t+)y 9) defined above. Then 

(1) (^(t_.t_|_), 3) is maximal globally hyperbolic. 

(2) There exists an uncountable number of analytic, vacuum, simply connected, 
non-equivalent extensions of {^^t_ t^j, g). 



44 PIOTR T. CHRUSCIEL, GREGORY J. GALLOWAY, AND DANIEL POLLACK 

6.2. Strong cosmic censorship. The strong cosmic censorship (SCC) problem 
concerns predictability: Indeed, a fundamental requirement of physically relevant 
equations is that solutions should be uniquely determined by initial data. So it is 
important to inquire about predictability in general relativity. 

In other words, we would like to know whether or not the solutions provided by 
Theorem 4.4 are unique. Now, it is easy to see that there can be no uniqueness 
unless some restrictions on the development are imposed: consider for example 
(-00, 1) X E", K X and (R x M") \ {(1, 0)} equipped with the obvious flat metric. 
All three space-times contain the spacelike surface ({0} x K",(5, 0), where 6 is the 
Euclidean metric on R". The first two are globally hyperbolic developments of 
the given initial data, but the third is not, as it is not globally hyperbolic. And 
obviously these are not isometric: e.g,. the second is geodesically complete, while 
the other two are not. So to guarantee uniqueness some further conditions are 
needed. 

The key existence and uniqueness theorem in this context is due to Choquet- 
Bruhat and Geroch [93] (compare [100,115,204]). Some terminology is needed: a 
space-time g) is said to be a development of an initial data set (Af , /i, K) if 
there exists an embedding i : M ^ such that «(Af) is a Cauchy surface for 
g), with 

i* g = h , 

and with K being the pull-back to M of the extrinsic curvature tensor (second 
fundamental form) of i{M). We will say that a development {^,g) is maximal 
globally hyperbolic if the following implication holds: if ip : ^ — >■ is an isomet- 
ric embedding of into (^', g'), and if V'(^) 7^ then is not globally 
hyperbolic. 

Note that wc are not imposing any field equations on (./#', g'). One could simi- 
larly define a notion of maximality within the class of vacuum space-times, but this 
would lead to a weaker statement of the Choquct-Bruhat - Geroch theorem, which 
for simplicity is presented in the smooth case: 

Theorem 6.2 (Existence of maximal globally hyperbolic developments [93]). For 
any smooth vacuum initial data (M, h, K) there exists a unique, up to isometric 
diffeomorphism, vacuum development (^,g), which is inextendible in the class of 
smooth globally hyperbolic Lorentzian manifolds. 

This theorem can be thought of as the equivalent of the usual ODE theorem of 
existence of maximal solutions. The generalization is, however, highly non-trivial 
because while the proof for ODEs deals with subsets of R, Theorem 6.2 deals with 
manifolds which are dynamically obtained by patching together local solutions. The 
main difficulty is to prove that the patching leads to a Hausdorff topological space. 
The argument makes use of Lorentzian causality theory, which in turn relies heavily 
on differentiability of the metric. To obtain a version of Theorem 6.2 with lower 
differentiability, as in Theorem 4.4 or in Remark 4.2, one would need to show that 
the relevant parts of causal theory can be repeated in the wider setting. 

While Theorem 6.2 is highly satisfactory, it does not quite prove what one wants, 
because uniqueness is claimed in the globally hyperbolic class only. But we have seen 
in Theorem 6.1 that there exist vacuum space-times with non-unique extensions of 
a maximal globally hyperbolic region. In such examples the space-time (^, g) of 
Theorem 6.2 is unique in the class of globally hyperbolic space-times, but it can be 
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extended in more than one way to strictly larger vacuum solutions. In such cases the 
extension always takes places across a Cauchy horizon, as defined in Section 2.3.3 
So one cannot expect uniqueness in general. However, it has been suggested by 
Penrose [319] that non-uniqueness happens only in very special circumstances. The 
following result of Isenberg and Moncricf [229,233,295] (compare [213]) indicates 
that this might indeed be the case: 

Theorem 6.3. Let {^,g) be a vacuum analytic space-time containing an analytic 
compact Cauchy horizon Jif. If the null geodesies threading Jff are closed, then the 
Cauchy horizon is a Killing horizon; in particular the isometry group of , g) is 
at least one-dimensional. 

The hypotheses of analyticity, compactness, and closed generators are of course 
highly restrictive. In any case it is conceivable that some kind of local isometrics 
need to occur in space-times with Cauchy horizons when those conditions are not 
imposed; indeed, all known examples have this property. But of course existence 
of local isometrics is a highly non-generic property, even when vacuum equations 
arc imposed [37], so a version of Theorem 6.3 without those undesirable hypotheses 
would indeed establish SCC. 

Whether or not Cauchy horizons require Killing vector fields, a loose mathe- 
matical formulation of strong cosmic censorship, as formulated in [114] following 
Moncricf and Eardley [294] and Penrose [319], is the following: 

Consider the collection of initial data for, say, vacuum or electro -vacuum 
space-times, with the initial data surface M being compact, or with 
asymptotically flat initial data (M, h, K) . For generic such data the 
maximal globally hyperbolic development is inextendible. 

Because of the difficulty of the strong cosmic censorship problem, a full under- 
standing of the issues which arise in this context seems to be completely out of 
reach at this stage. There is therefore some interest in trying to understand that 
question under various restrictive hypotheses, e.g., symmetry. The simplest case, of 
spatially homogeneous space-times, has turned out to be surprisingly difficult, be- 
cause of the intricacies of the dynamics of some of the Bianchi models discussed in 
Section 6.7, and has been settled in the affirmative in [143] (compare Theorem 6.9 
below). 

6.2.1. Gowdy toroidal metrics. The next simplest case is that of Gowdy metrics on 
T3 := X X S^: by definition, 

(6.4) g = e'-^-^^/'^-er'^^dr^ +de'^) + e-^ePda^ + 2ePQdcrdS+{ePQ^ + e-P)dS% 

where r e R and (6,(7,5) are coordinates on T'^, with the functions P,Q and 
A depending only on r and 9. The metric of a maximal globally hyperbolic 
U{1) X J7(l)-symmetric vacuum space-time with T^^-Cauchy surfaces can be glob- 
ally written [113] in the form (6.4) provided that the twist constants vanish: 

(6.5) ca := e^p^sX-^X^V^Xi = , a = 1, 2 , 

where the Xa's are the Killing vectors generating the f/(l) x U{1) action. The 
condition ci = C2 = is equivalent to the requirement that the family of planes 
span{Xi,X2}^ is integrable. 
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For metrics of the form (6.4), the Einstein vacuum equations become a set of 
wave-map equations 



which is a consequence of (6.9) and of periodicity in 9. The metric function A is 
obtained by integrating (6.8)-(6.9). Global existence of solutions to (6.6)-(6.7) was 
proved in [292] when the initial data are given on a hypersurface {r = const}, and 
in [113] for general U{1) x [/(l)-symmetric Cauchy surfaces. 

The question of SCC in this class of metrics has been settled by Ringstrom, who 
proved that the set of smooth initial data for Gowdy models on that do not 
lead to the formation of Cauchy horizons contains a set which is open and dense 
within the set of all smooth initial data. More precisely, Ringstrom's main result 
(see [334, 335] and references therein) is the following: 

Theorem 6.4. Let tq G M and let S = {((5(to), P(to), Qr (tq), ^t(7'o))} be the set 
of smooth initial data for (6.6)-(6.7) satisfying (6.10). There is a subset Q of S 
which is open with respect to the x topology, and dense with respect to the 
C°° topology, such that the space-times of the form (6.4) corresponding to initial 
data in Q are causally geodesically complete in one time direction, incomplete in the 
other time direction, and the Kretschmann scalar, Rap^sR"^'*^ , becomes unbounded 
in the incomplete direction of causal geodesies. 

This result does indeed establish SCC in this class of metrics: to see that the 
resulting space-times are inextendible in the category of manifolds with 
Lorentzian metrics, note that the existence of any such extension would imply 
existence of geodesies which are incomplete in the original space-time, and along 
which every curvature scalar is bounded. 

Theorem 6.4 is complemented by the results in [84,134,293], where infinite di- 
mensional families of (nongcneric) solutions which are extendible across a Cauchy 
horizon are constructed. 

The key to the understanding of the global structure of the Gowdy space-times 
is the analysis of the behavior of the functions P and Q as r — > ±oo. The asymp- 
totic behavior of those functions, established by Ringstrom, can then be translated 
into statements about the behavior of the space-time geometry as those limits are 
approached. A central element of the proof is the existence of a velocity function 



Essential steps in Ringstrom's analysis are provided by the work on Fuchsian PDEs 
of Kichenassamy and Rendall [241,330], as well as the study of the action of Geroch 




Initial data on for P and Q have to satisfy an integral constaint. 



(6.10) 
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transformations by Rendall and Weaver [331] (compare [84]). See also [137] for the 
related problem of an exhaustive description of Cauchy horizons in those models. 

6.2.2. Other U (1) x U {!)- symmetric models. The existence of two Killing vectors is 
also compatible with S^, L{p,q) ("lens" spaces), and x topologies. Thus, to 
achieve a complete understanding of the set of spatially compact initial data with 
precisely two Killing vectors one needs to extend Ringstrom's analysis to those cases. 
There is an additional difficulty that arises because of the occurrence of axes of 
symmetry, where the (l + l)-reduced equations have the usual singularity associated 
with polar coordinates. Nevertheless, in view of the analysis by Christodoulou and 
Tahvildar-Zadeh [109,110] (see also [113]), the global geometry of generic maximal 
globally hyperbolic solutions with those topologies is reasonably well understood. 
This leads one to expect that one should be able to achieve a proof of SCC in those 
models using simple abstract arguments, but this remains to be seen. 

Recall, finally, that general models with two Killing vectors Xi and X2 on 
have non-vanishing twist constants (6.5). The Gowdy metrics are actually "zero 
measure" in the set of all U{1) x U{1) symmetric metrics on T'^ because Ca = for 
the Gowdy models. The equations for the resulting metrics are considerably more 
complicated when the Ca's do not vanish, and only scant rigorous information is 
available on the global properties of the associated solutions [40,235,327]. It seems 
urgent to study the dynamics of those models, as they are expected to display [41] 
"oscillatory behavior" as the singularity is approached, in the sense of Section 6.8. 
Thus, they should provide the simplest model in which to study this behavior. 

6.2.3. Spherical symmetry. One could think that the simplest possible asymptoti- 
cally flat model for studying the dynamics of the gravitational field will be obtained 
by requiring spherical symmetry, since then the equations should reduce to wave 
equations in only two variables, t and r. Unfortunately, for vacuum space-times this 
turns out to be useless for this purpose because of Birkhoff's theorem [49], which 
asserts that spherically symmetric vacuum metrics are static. So, if one wishes to 
maintain spherical symmetry, supplementary fields are needed. The case of a scalar 
field was studied in a series of intricate papers over 13 years by Christodoulou, begin- 
ning with [102], and culminating in [104] with the verification of the strong cosmic 
censorship conjecture within the model. Christodoulou further established "weak 
cosmic censorship" in this class, an issue to which we return in the next section, 
and exhibited non-generic examples for which the conclusions of these conjectures 
fail [103]. 

The situation changes when electromagnetic fields are introduced. The analysis 
by Dafermos [151,152] of the spherically symmetric Einstein-Maxwell-scalar field 
equations yields a detailed picture of the interior of the black hole for this model, 
in terms of initial data specified on the event horizon and on an ingoing null hy- 
persurface. When combined with the work by Dafermos and Rodnianski [153] on 
Price's law, one obtains the following global picture: initial data with a compactly 
supported scalar field, and containing a trapped surface (see Section 7.2 below), 
lead to space-times which either contain a degenerate (extremal) black hole, or de- 
velop a Cauchy horizon, with a space-time metric that can be continued past this 
horizon in a , hut not manner. It seems that not much is known about the 
properties of the degenerate solutions, which are presumably non-generic; it would 
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be of interest to clarify tliat. In any case, the work shows that strong cosmic cen- 
sorship holds within the class of nondegenerate solutions with trapped surfaces, at 
the level, leaving behind the perplexing possibility of continuous extendability 
of the metric. 

The reader is referred to [3,114,329] and references therein for further reading 
on sec. 

6.3. Weak cosmic censorship. The strong cosmic censorship conjecture is an 
attempt to salvage predictability of Einstein's theory of gravitation. There exists 
a variant thereof which addresses the fact that we do not seem to observe any of 
the singularities that are believed to accompany gravitational collapse. The hope is 
then that, generically, in asymptotically flat space-times, any singular behavior that 
might form as a result of gravitational collapse, such as causality violations, lack of 
predictability, or curvature singularities, will be clothed by an event horizon. For 
this, one introduces the notion of future null infinity, which is an idealized boundary 
attached to space-time that represents, loosely speaking, the end points of null 
geodesies escaping to infinity. (In stationary situations this is closely related to the 
region ^cxt of (3.27).) The black hole event horizon is then the boundary of the past 
of null infinity; compare (3.28) and (3.29). One then wishes the part of the space- 
time that lies outside the black hole region to be well-behaved and "sufficiently 
large" . This is the content of the weak cosmic censorship conjecture, originally due 
to Penrose [319], as made precise by Christodoulou [105]: for generic asymptotically 
flat initial data, the maximal globally hyperbolic development has a complete future 
null infinity. Heuristically this means that, disregarding exceptional sets of initial 
data, no singularities are observed at large distances, even when the observations are 
continued indefinitely. One should remark that, despite the names, the strong and 
weak cosmic censorship conjectures are logically independent; neither follows from 
the other. Note also that some predictability of Einstein's theory would be salvaged 
if strong cosmic censorship failed with weak cosmic censorship being verified, since 
then the failure of predictability would be invisible to outside observers. 

Both cosmic censorship conjectures arc intimately related to the issue of grav- 
itational collapse, the dynamical formation of black holes and singularities, first 
observed for a homogeneous dust model by Oppenheimer and Snyder in 1939 [309], 
visualized in Figure 6.1. 

So far the only complete analysis of weak cosmic censorship in a field theo- 
retical model is that of the spherically symmetric scalar field model studied by 
Christodoulou [103,104], already mentioned in Section 6.2.3. 

6.4. Stability of vacuum cosmological models. Not being able to understand 
the dynamics of all solutions, one can ask whether some features of certain par- 
ticularly important solutions persist under small perturbations of initial data. For 
example, will geodesic completeness still hold for space-times arising from small 
perturbations of Minkowskian initial data? Or, will a global, all encompassing, 
singularity persist under perturbations of Bianchi IX initial data (see Section 6.7). 
Such questions are the object of stability studies. 

6.4.1. U{1) symmetry. Our understanding of models with exactly one Killing vector 
is dramatically poorer than that of U{1) x U{1) symmetric space-times. Here one 
only has stability results, for small perturbations within the U{1) isometry class 
in the expanding direction ("away from the singularity"): In [87] Choquet-Bruhat 
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considers U{1) symmetric initial data {h, K) for the vacuum Einstein equations on 
a manifold of the form M x S^, where M is a compact surface of genus g > I. It 
is assumed that trft,A' is constant, and that {h, K) are sufficiently close to (ft-o, K^)), 
where /iq is a product metric 

^0 = 7 + dx^ , 

with 7 being a metric of constant Gauss curvature on M , and with Kq proportional 
to ho- The sign of the trace of Kq determines an expanding time direction and a 
contracting one. Under those conditions, Choquet-Bruhat proves that the solution 
exists for an infinite proper time in the expanding direction. The analysis builds 
upon previous work by Choquet-Bruhat and Moncrief [99] , where a supplementary 
polarization condition has been imposed. Not much is known in the contracting 
direction in the C/(l)-symmetric models (see, however, [232]), where "mixmaster 
behavior" is expected [39,43]; compare [44]. 

6.4.2. Future stability of hyperbolic models. The proof of the above result bears 
some similarity to the future stability theorem of Andersson and Moncrief [13], as 
generalized in [11], for spatially compact hyperbolic models without any symmetries. 
Those authors consider initial data near a negatively curved compact space form, 
with the extrinsic curvature being close to a multiple of the metric, obtaining future 
geodesic completeness in the expanding direction. The control of the solution is 
obtained by studying the Bel-Robinson tensor and its higher-derivatives analogues. 
A striking ingredient of the proof is an elliptic-hyperbolic system of equations, used 
to obtain local existence in time [12]. 

6.5. Stability of Minkowski space-time. 



"'^^See the discussion after Theorem 6.8, and Section 6.8. 
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6.5.1. The Christodoulou-Klainerman proof. One of the flagship results in mathe- 
matical general relativity is nonlinear stability of Minkowski space-time, first proved 
by Christodoulou and Klaincrman [107]. One starts with an asymptotically flat vac- 
uum initial data set {h,K) on M.^. Under standard asymptotic flatness conditions, 
for {h,K) sufficiently close to Minkowskian data, the maximal globally hyperbolic 
development {^^,g) of the data contains a maximal hypcrsurface, i.e., a hypersur- 
face satisfying lihK = 0; this follows from the results in [22,28,108]. So without 
loss of generality one can, in the small data context, assume that the initial data 
set is maximal. 

The precise notion of smallness needed for the Christodoulou-Klainerman theo- 
rem is defined as follows: For p G S ~ M"^, a > 0, consider the quantity 

» 1 2 

(6.11) 0(a,p) = / {^(4 + a2)^+i|V^Ric|2 + ^(^ + a')'\^' K\^}d^lg , 

•^^ e=o e=i 

where dp is the geodesic distance function from p, Ric is the Ricci tensor of the 
metric g, dfig is the Riemannian measure of the metric g and V is the Riemannian 
connection of g. Let 

Q» = inf Q{a,p) . 

Christodoulou and Klaincrman prove causal geodesic completeness of (^, g) pro- 
vided that is sufficiently small. The proof proceeds via an extremely involved 
bootstrap argument involving a foliation by maximal hypersurfaces St together 
with an analysis of the properties of an optical function u. In the context here this 
is a solution of the eikonal equation 

g'^^daudpu = , 

the level sets C„ of which intersect St in spheres which expand as t increases. We 
have: 

Theorem 6.5 (Global Stability of Minkowski space-time). Assume that (M, h, K) 
is maximal, with^^ 

(6.12) /ly = 6,, + 03{r-'/^), lU, ^ 02(r^'/') . 

There is an e > such that if < e, then the maximal globally hyperbolic devel- 
opment {^,g) of {AI, h, K) is geodesically complete. 

The above version of Theorem 6.5 is due to Bieri [48]. The original formulation 
in [107] assumes moreover that 

(6.13) h^ {l + 2]\I/r)S + Oiir-^/^), K = osir'^/^) , 

and in the definition (6.11) a term involving K with £ = is added. 

By definition, asymptotically fiat initial data sets approach the Minkowskian 
ones as one recedes to infinity. One therefore expects that at sufficiently large dis- 
tances one should obtain "global existence" , in the sense that the maximal globally 
hyperbolic development contains complete outgoing null geodesies. This question 
has been addressed by Klainerman and Nicolo [244-246]; the reader is referred to 
those references for precise statements of the hypotheses made: 



-•^^A function / on M is oj.(r if r^+'V'/ — 5> as r — >■ oo for all i = 0, . . . fc. 
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Theorem 6.6. Consider an asymptotically flat initial data set {M,h, K), with 
maximal globally hyperbolic development (^,g). Let fir denote a conditionally 
compact domain bounded by a coordinate sphere Sr C ^cxt- There exists R > 
such that for all r > R the generators of the boundary dJ'^{Qr) of the domain of 
influence J~^{Qr) of fir are future- complete. 

Both m [107] and in [245] one can find detailed information concerning the 
behavior of null hypcrsurfaces as well as the rate at which various components 
of the Ricmann curvature tensor approach zero along timelike and null geodesies. 

6.5.2. The Lindblad-Rodnianski proof. A completely new proof of stability of Minkowski 
space-time has been given by Lindblad and Rodnianski [267, 268] . The method 
provides less detailed asymptotic information than [107] and [245] on various quan- 
tities of interest but is much simpler. The argument is flexible enough to allow the 
inclusion of a scalar field, or of a Maxwell field [273, 274] (compare [48] for an anal- 
ysis along the lines of the Christodoulou-Klainerman approach) , and generalizes to 
higher dimensions [90]. Further it allows the following, less restrictive than that 
in [107,245], asymptotic behavior of the initial data, for some a > 0: 

(6.14) h= {l + 2m/r)d + 0{r-^-°') , K ^ 0{r-^~°') . 

Lindblad and Rodnianski consider the Einstein-Maxwell equations with a neutral 
scalar field: 

R 

2' 

with 



(6.15) Rij,y - —9^,1, = Tfj^^ + Tfj,y , 



T, 



The initial data are prescribed on M", so that the Maxwell field F has a global 
potential A, F^^ ~ d^A^ — d^Ai^^. The matter field equations read 

(6.16) D^F^''' = , D,,^. = . 

The initial data, denoted by (h, K, A, E, -0o, "^i) (where, roughly speaking, A is the 
initial value for the Maxwell potential and E is the initial value for the electric 
field), satisfy the following asymptotic conditions, for r = [x] — > oo, with some 
a > 0: 



hij 



3.17) 



(l + 2f)%+0(r-i-") , forn = 3, 
+ 0{r^-°') , for n > 4 , 



A = o{r^-°') , kij = o{r-'^-°') , i; = 0(r-T^-") , 



One also supposes that the Einstein-Maxwell constraint equations hold initially: 

r R = k'^k.j - k^kji + 2E,k + f,,f'^ + jWol' + l^il' , 

(6.18) <^ V^k,j-V,kj^ = FojF^ +V,,Po,Pi , 

where R is the scalar curvature of the metric h. The strategy is to impose globally 
the wave coordinates condition 

(6.19) 5^(ff''V|det<?l) =0 Vi. = 0,...,n, 
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as well as the Lorenz gauge for the electromagnetic potential Af^,, 



(6.20) 5,, (^Vldetffl^^j =0. 
Letting Dg = g^'^d^d^, the dynamical equations take the form 

(6.21) nj A, I = I ] \ ^ 

where the source terms S^u and are bilinear in the derivatives of the fields, with 
coefiicients depending upon the metric. The initial data are decomposed as 

\{r/t)x{ry-f5^.u forn = 3 , 



3.22) /i^, = V-C' ^ithCW 



for n > 4 , 



where x G C"^ is any function such that x{s) equals 1 for s > 3/4 and for s < 1/2. 
The proof relies heavily on the structure of the nonlinear terms in wave coordinates. 

Recall that there exists an extensive literature on wave equations in 3 + 1 di- 
mensions with nonlinearities satisfying the null condition [242,243], but the nonlin- 
earities that arise do not satisfy that condition. The argument works only because 
different components of h can be treated on a different footing. Indeed, for solu- 
tions of the wave equation on Minkowski space-time, the derivatives in directions 
tangent to the light cones decay faster than the transverse ones. But the wave 
coordinates condition (6.19) can be used to express the transverse derivatives of 
some components of g^^, in terms of tangential derivatives of the remaining ones. 
This provides control of the nonlinearities. 

We also note the small data global existence results of [220,262] on M""'"^, n > 4,, 
and of [101] for odd n > 5. The structure conditions there are general enough to 
cover the Einstein equations in wave coordinates, but the assumptions on the fall-off 
of initial data exclude non-trivial solutions of the vacuum constraint equations^^ 

We have: 

Theorem 6.7. Consider smooth initial data {h,K,A,E,^Q,^i) on R", n > 3, 
satisfying (6.17) together with the Einstein- Maxwell constraint equations. Let N G 
N, suppose that := N + [^^±2] _ 2 > 6 2^^], and set 

(6.23) 

En.AO)- E (ll(l + r)V2+7+mvvX||i.-f ||(H-r)V2+7+mv^X||i, 

0<\I\<N„ 

+ 11(1 + r)i/2+^+I^IVV^i||i. + 11(1 + r)i/2+^+I^IV^^||2, 

+ 11(1 + ^)1/2+7+|/|vvVo||l^ + 11(1 + r)i/2+7+|/|vVi|U.) . 

Let m be the ADM mass of h. For every 70 > there exists Eq > 0, with 70 (£0) ~^ 
as Sq — > 0, such that if 



(6.24) ^/eZMO) + m < eo, 

^^In [220, 262] compactly supported data arc considered. In the theorem for general quasi- 
hnear systems given in [101] the initial data are in a Sobolev space which requires fall-off at 
infinity faster than j--"— 3/2 jjj both cases the positive energy theorem implies that such initial 
data lead to Minkowski space-time. 
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for some 7 > 70, then the maximal globally hyperbolic development of the initial 
data is geodesically complete. 

The proof by Lindblad and Rodnianski is an ingenious and intricate analysis 
of the couphng between the wave-coordinates gauge and the evolution equations. 
One makes a clever guess of how the fields decay in space and time, encoded in the 
following weighted energy functional, 

(6.25) 

'■(i)^ sup y / {\dZ'h'\^ + \dZ'A\^ + \dZ'i;f)w{q)d''x , 
o<T<t r^^J—f.^^, Jt 



Matter / 



' Z<£2r,\I\<N„ 



where 

(6.26) w{q) 



1 + (1 + M)'+'^ q>0 

l + (l + ky|)-'^ '7<o, 



with q~r~t,iJL>{) and < 7 < 1 . Here denotes the collection of the 
following generators of the conformal Lorentz group, first used to study the decay 
of solutions of the Minkowski wave equation by Klainerman [242] : 

(6-27) 9., 

One argues by continuity: one chooses < (S < -j , and one considers the maximal 
time T so that the inequality 

(6.28) •^^j^„"""'(t) < 2Cw„e'(l + tf^ 

holds for Q < t < T . A sophisticated method, using the Klaincrman-Sobolev 
inequalities [242], together with a new weighted energy inequality, allows one to 
show that (6.28) holds for < i < 2^ with a smaller constant on the right-hand- 
side, contradicting maximality of T , and thus proving global existence. 

A long standing question in the study of asymptotically flat space-times is that 
of the existence of an asymptotic expansion of the metric as one recedes to infinity 
along outgoing null cones, see [180,182,318]. Neither the analysis of [267,268], nor 
that in [107,245], provides sufficient information. It would be of interest to clarify 
that. 



6.6. Towards stability of Kerr: wave equations on black hole backgrounds. 

Since the pioneering work of Christodoulou and Klainerman on stability of Minkowski 
space-time, many researchers have been looking into ways to address the question 
of stability of Kerr black holes. The first naive guess would be to study stability 
of Schwarzschild black holes, but those cannot be stable since a generic small per- 
turbation will introduce angular momentum. The current strategy is to study, as a 
first step, linear wave equations on black hole backgrounds, with the hope that suf- 
ficiently robust linear decay estimates can be bootstrapped to produce a nonlinear 
stability proof. Due to limited space we will not review those results, referring the 
reader to recent important papers on the subject [4,52,53,154-156,163,278,361,362], 
see also [176] and references therein. 
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Table 6.1. Lie groups of Bianchi class A. 



Bianchi type 


ni 


?^2 




Simply connected group 




I 











Abelian 




II 


+ 








Heisenberg 




VIo 





+ 




Sol (isometries of the Minkowski plane 1 




VIIo 





+ 


+ 


universal cover of Euclid (isometries of I 




VIII 




+ 


+ 


universal cover of SL{2,R) 




IX 


+ 


+ 


+ 


SU{2) 





6.7. Bianchi A metrics. Another important example of the intricate dynamical 
behavior of solutions of the Einstein equations is provided by the "Bianchi A" 
vacuum metrics. The key insight provided by these space-times is the supposedly 
chaotic behavior of large families of metrics in this class when a singularity is 
approached. This dynamics has been conjectured to be generic; we will return 
to this issue in Section 6.8. As will be seen shortly, in Bianchi A space-times 
the Einstein evolution equations reduce to a polynomial dynamical system on an 
algebraic four-dimensional submanifold of K^. The spatial parts of the Bianchi 
geometries provide a realization of six, out of eight, homogeneous geometries in 
three dimensions which form the basis of Thurston's geometrization program. 

For our purposes here we define the Bianchi space-times as maximal globally 
hyperbolic vacuum developments of initial data which are invariant under a simply 
transitive group of isometries. Here the transitivity of the isometry group is meant 
at the level of initial data, and not for the space-time. The name is a tribute 
to Bianchi, who gave the classification of three dimensional Lie algebras which 
underline the geometry here. These metrics split into two classes, Bianchi A and 
Bianchi B, as follows: Let G be a 3-dimcnsional Lie group, and let Zi, i = 1,2,3 
denote a basis of left-invariant vector fields on G. Define the structure constants 
lij by the formula 

The Lie algebra and Lie group are said to be of class A if ^ik' = 0; class B are the 
remaining ones. The classes A and B correspond in mathematical terminology to 
the unimodular and non-unimodular Lie algebras. A convenient parameterization 
of the structure constants is provided by the symmetric matrix n*^ defined as 

(6.29) n'J" = ^jki^'e^^"'. 

This implies jij'' ~ ^ijmn'^™- The Bianchi A metrics are then divided into six 
classes, according to the eigenvalues of the matrix n'-', as described in Table 6.1. 
For the Bianchi IX metrics, of particular interest to us here, the group G is SU (2). 
Thus, the Taub metrics discussed in Section 6.1 are members of the Bianchi IX 
family, distinguished by the existence of a further t/(l) factor in the isometry group. 

Let G be any three-dimensional Lie group, the Lie algebra of which belongs 
to the Bianchi A class. (The G's are closely related to the Thurston geometries, 
see Table 6.1; compare [3, Table 2]). Denote by {cr'} the basis dual to {Zi}. It is 
not too difficult to show that both A and B Bianchi metrics can be globally written 
as 



(6.30) 



g = -dt'^ + hij {t)a'a^ , t e I , 
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with a maximal time interval /. 

There are various ways to write the Einstein equations for a metric of the form 
(6.30). We use the formalism introduced by Wainwright and Hsu [368], which has 
proven to be most useful for analytical purposes [328,332,333], and we follow the 
presentation in [332] . Let 

(Trj = A% - ifr/iA'/iy , 9 := tihK , 

be the trace-free part of the extrinsic curvature tensor of the level sets of t. Away 
from the (isolated) points at which 9 vanishes, one can introduce 

Nij = nij/9 , 



Sij — Bij -rB\Sij 



Set S+ = 1(^22 + S33) and E_ = ^3(1122 - ^33)72. If we let Ni be the eigenvalues 
of Nij , the vacuum Einstein equations (a detailed derivation of which can be found 
in [332]) lead to the following autonomous, polynomial dynamical system 

N[ = (g-4E+)iVi, 

= (g + 2E+ + 2V3I]_)Af2 , 

(6.31) = (g + 2E+ - 2V3I]_)A^3 , 

K = -(2-g)S+-35+ , 

E'_ = -(2-g)I]_ -3S'_ , 

where a prime denotes derivation with respect to a new time coordinate r defined 

by 



(6.32) 
Further, 



dt 
d^ 



q = 2(E; + S]i) 



(6.33) S+ = ^[{N2 - Nsf ^ Ni{2Ni - N2 - N3)] , 

S- = ^{N3~N2)iN,-N2-N3) . 
The vacuum constraint equations reduce to one equation, 

(6.34) E^ + E^ + ^ [N^ + iV| + - 2{NiN2 + iV2iV3 + N^Ni)] = 1 . 

The points (iVi, A^2, -^3, E_) can be classified according to the values of iVi, A^2, -^3 
in the same way as the n"s in Table 6.1. The sets iV^ > 0, A^i < and Ni = are in- 
variant under the flow determined by (6.31), and one can therefore classify solutions 
to (6.31)-(6.34) accordingly. Bianchi IX solutions correspond, up to symmetries of 
the system, to points with all TV^'s positive, while for Bianchi VIII solutions one 
can assume that two Ni 's are positive and the third is negative. 
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Points with Ni ~ N2 ^ N3 ^ correspond to Bianchi / models. The associated 
vacuum metrics were first derived by Kasner, and take the form 

3 

(6.35) ds^ = -df + ^ t'^P'dx^ ® dx\ pi + p2 + ps = pi + pj + pi = I ■ 
1=1 

An important role in the analysis of (6.31) is played by the Kasner circle, defined 
as the set {q = 2}. These points belong to the configuration space, as determined 
by (6.34), for Bianchi I models, but the equation g = 2 is incompatible with (6.34) 
for Bianchi IX metrics. Nevertheless, we shall see shortly that the Kasner circle 
plays an essential role in the analysis of the Bianchi IX dynamics. 

The set E_ = 0, N2 = N3, together with its permutations, is invariant under 
the flow of (6.31)-(6.34). In the Bianchi IX case these are the Taub solutions. In 
the Bianchi VIII case the corresponding explicit solutions, known as the NUT 
metrics, have been found by Newman, Tamburino and Unti [301], and they exhibit 
properties similar to the Bianchi IX Taub solutions discussed in Section 6.1. 

The Lo-limit of an orbit 7 of a dynamical system is defined as the set of accumu- 
lation points of that orbit. In [332,333], Ringstrom proves the following: 

Theorem 6.8. The ui-limit set of each non-NUT Bianchi VIII orbit contains at 
least two distinct points on the Kasner circle. Similarly, non-Taub-NUT Bianchi 
IX orbits have at least three distinct uj-limit points on the Kasner circle. 

The picture which emerges from a numerical analysis of (6.31) (see [42, 147] 
and references therein) is the following: Every non-Taub-NUT Bianchi IX orbit 
approaches some point on the Kasner circle; there it performs a "bounce" , after 
which it eventually approaches another point on the Kasner circle, and so on. 
Theorem 6.8 establishes the validity of this picture. The numerical analysis further 
suggests that generic orbits will have a dense w-limit set on the Kasner circle; this 
is compatible with, but does not follow from, Ringstrom's analysis. It has been 
argued that the map which associates to each bounce the nearest point on the 
Kasner circle possesses chaotic features; this is at the origin of the "mixmaster 
behavior" terminology, sometimes used in this context. Major progress concerning 
this issue has been achieved in [32,266,325], where existence of orbits exhibiting 
the above behaviour has been established, but the question of what happens for 
all, or for generic orbits remains open. 

The following result of Ringstrom [332] provides further insight into the geometry 
of Bianchi IX space-times: 

Theorem 6.9. In all maximal globally hyperbolic developments {^,g) of non- 
Taub-NUT Bianchi IX vacuum initial data or of non-NUT Bianchi VIII vacuum 
initial data the Kretschmann scalar 

is unbounded along inextendible causal geodesies. 

Note that the observation of curvature blow-up provides a proof, alternative to 
that of [143], of the non-existence of Cauchy horizons in generic Bianchi IX models. 

We close this section by mentioning that no similar rigorous results are known 
concerning the global dynamical properties of Bianchi models of class B; com- 
pare [210]. 
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Figure 6.2. A few "bounces" in a typical Bianchi IX orbit; figures 
and numerics by Woei-Chet Lim. The vertical axis represents A^i 
(red), N2 (green), 7V3 (blue), with only the biggest of the iV^'s 
plotted. The Kasner circle and the triangle for the Kasner billiard 
in the (E+, S_)-plane are shown. The projected trajectories can 
be seen to approach the billiard ones. 

6.8. The mixmaster conjecture. The most important question in the study of 
the Cauchy problem is that of the global properties of the resulting space-times. 
So far we have seen examples of geodesically complete solutions (e.g., small per- 
turbations of Minkowski space-time), or all-encompassing singularities (e.g., generic 
Bianchi models), or of Cauchy horizons (e.g., Taub-NUT metrics). The geodesically 
complete solutions are satisfying but dynamically uninteresting, while the strong 
cosmic censorship conjecture expresses the hope that Cauchy horizons will almost 
never occur. So it appears essential to have a good understanding of the remaining 
cases, presumably corresponding to singularities. Belinski, Khalatnikov and Lifs- 
chitz [39] suggested that, near singularities, at each space point the dynamics of the 
gravitational field resembles that of generic Bianchi metrics, as described in Sec- 
tion 6.7. Whether or not this is true, and in which sense, remains to be seen; in any 
case the idea, known as the BKL conjecture, provided guidance — and still does — 
to a significant body of research on general relativistic singularities; see [21,159,165] 
and references therein. This then leads to the mathematical challenge of making 
sense of the associated slogan, namely that the singularity in generic gravitational 
collapse is spacelike, local, and oscillatory. Here spacelike is supposed to mean that 
strong cosmic censorship holds. The term local refers to the idea that, near generic 
singularities, there should exist coordinate systems in which the metric asymptotes 
to a solution of equations in which spatial derivatives of appropriately chosen fields 
have been neglected. Finally, oscillatory is supposed to convey the idea that the 
approximate solutions will actually be provided by the Bianchi IX metrics. 



°The resulting truncated equations should then presumably resemble the equations satisfied 
by spatially homogeneous metrics. However, different choices of quantities which are expected 
to be time-independent will lead to different choices of the associated notion of homogeneity; for 
instance, in [39] the types Bianchi VIII and IX are singled out; the notion of genericity of those 
types within the Bianchi A class is read from Table 6.1 as follows: "something that can be non- 
zero is more generic than something that is" . On the other hand, the analysis in [206] seems to 
single out Bianchi VI_i/g metrics. 
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The main rigorous evidence for a relatively large class of vacuum^^ space-times 
with singularities which are spacelike and local in the sense described above is 
Ringstrom's Theorem 6.4, describing generic Gowdy metrics, but the resulting sin- 
gularities are not oscillatory. This is not in contradiction with the conjecture, since 
the Gowdy metrics arc certainly not generic, whether in the space of all metrics, or 
in the space of U{1) x U{1) symmetric ones: As mentioned in Section 6.4.1, generic 
[/(I) X t7(l) metrics have non-vanishing twist constants as defined by (6.5). The 
numerical studies of [41] suggest that the switching-on of the twist constants will 
indeed generically lead to some kind of oscillatory behavior. 

In fact, BKL put emphasis on Bianchi IX models, while some other authors 
seem to favor Bianchi VI_i/g, or not-necessarily Bianchi, oscillations [39,43,159, 
206,366]. It has moreover been suggested that the oscillatory behavior disappears 
in space-time dimensions higher than ten [159, 161], and large families of non- 
oscillatory solutions with singularities have indeed been constructed in [160]. This 
leads naturally to the following, somewhat loose, conjecture: 

Conjecture 6.10 (Mixmaster conjecture). Let n + 1 < 10. There exist open sets 
of vacuum metrics for which some natural geometric variables undergo oscillations 
of increasing complexity along inextendible geodesies of unbounded curvature. 

The BKL conjecture would thus be a more precise version of the above, claiming 
moreover gencricity of the behavior, and pointing out to the Bianchi dynamics as 
the right model. Those properties are so speculative that we decided not to include 
them in Conjecture 6.10. 

The only examples so far of oscillatory singularities which are not spatially ho- 
mogeneous have been constructed by Berger and Moncrief [44]. There, a solution- 
generating transformation is applied to Bianchi IX metrics, resulting in non-ho- 
mogeneous solutions governed by the "oscillatory" functions arising from a non- 
Taub Bianchi IX metric. The resulting metrics have at least one but not more 
than two Killing vectors. The analysis complements the numerical evidence for 
oscillatory behavior in U{1) symmetric models presented in [45]. 

7. MARGINALLY TRAPPED SURFACES 

There have been some interesting recent developments at the interface of space- 
time geometry and the theory of black holes associated with the notion of marginally 
outer trapped surfaces. Let S be a co-dimension two spacelike submanifold of a 
space-time ^ . Under suitable orientation assumptions, there exist two families of 
future directed null geodesies issuing orthogonally from E. If one of the families has 
vanishing expansion along E, then E is called a marginally outer trapped surface 
(or an apparent horizon) . The notion of a marginally outer trapped surface was in- 
troduced early on in the development of the theory of black holes, as the occurrence 
of the former signals the presence of the latter. More recently, marginally outer 
trapped surfaces have played a fundamental role in quasi-local descriptions of black 
holes, and have been useful in numerical simulations of black hole space-times; see 
e.g. [17]. Marginally outer trapped surfaces arose in a more purely mathematical 
context in the work of Schoen and Yau [347] concerning the existence of solutions 
to the Jang equation, in connection with their proof of the positivity of mass. 



^^See, however, [14, 160] for a class of space-times with sources; [160] also covers vacuum in 
space dimensions n > 10. 
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Mathematically, marginally outer trapped surfaces may be viewed as space-time 
analogues of minimal surfaces in Riemannian manifolds. Despite the absence of 
a variational characterization like that for minimal surfaces, marginally trapped 
surfaces have recently been shown to satisfy a number of analogous properties, cf., 
in particular, [7-10,16,166,190] and the recent review [5]. The aim of this section 
is to describe some of these mathematical developments. 

7.1. Null hypersurfaces. Each family of null geodesies issuing orthogonally from 
S, as described above, forms a smooth null hypersurface near E. It would be 
useful at this stage to discuss some general aspects of such hypersurfaces. Null 
hypersurfaces have an interesting geometry, and play an important role in general 
relativity. In particular, as we have seen, they represent horizons of various sorts, 
such as the event horizons discussed in Section 3. 

Let (^"+^, p) be a space-time, with n > 2. A smooth null hypersurface in ^ is 
a smooth co-dimension one submanifold ^ = of ^ such that the restriction of 
g to each tangent space of -yV is degenerate. This, together with the Lorentz 

signature, implies that there is a unique direction of degeneracy in each tangent 
space . Thus, every null hypersurface .yV comes equipped with a smooth 

future directed null vector field K {g{K, K) ~ 0) defined on, and tangent to ^ 
such that the normal space of K at each p € yV coincides with the tangent space of 
^ at p, i.e., = TpjV for all p € . Tangent vectors to jV , transverse to K , 
arc then necessarily spacelike. The null vector field K associated to ,jV is unique 
up to positive pointwise rescaling. However, there is, in general, no canonical way 
to set the scaling. 

Two simple examples arise in Minkowski space M^'". The past and future cones 
dl^ (p) and 9/^ (p) are smooth null hypersurfaces away from the vertex p. Each 
nonzero null vector v £ TpM}'" determines a null hyperplane H = {q £ R^'" : 
r]{pq, v) ~ 0}, where 77 is the Minkowski metric, and pq is the tangent vector at p 
representing the displacement from p to q. 

It is a fundamental fact that the integral curves of K are null geodesies, though 
perhaps not affincly parameterized - this will depend on the scaling of K . Thus ^ 
is ruled by null geodesies, called the null generators of . For example the future 
cone jV = dl~^{p) \ {p} in Minkowski space is ruled by future directed null rays 
emanating from p. 

The null expansion scalar 9 of ^ with respect to if is a smooth function on yi^ 
that gives a measure of the average expansion of the null generators of y^ towards 
the future. In essence, 9 is defined as the divergence of the vector field K along y^. 
To be precise, given p £ yy, let n„_i be a co-dimension one subspace of TpyV"^ 
transverse to Kp. The metric restricted to n„_i, will be positive definite. Let 
{ei, 62, e„_i} be an orthonormal basis for n„_i with respect to g. Then 9 at p 
is defined as. 



Interestingly, due to the fact that K is null, this value is independent of the choice 
of transverse subspace n„_i, as well as of the choice of an orthonormal basis for 
n„_i; and so the expansion scalar 9 is well defined. 



n-l 



(7.1) 




^''There seems to be no analogue of the area functionaL 
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While 6 depends on the choice of K, it does so in a simple way. As easily 
follows from Equation (7.1), a positive rescaling of K rescales 9 in the same way: 
\i K = fK then = JO. Thus the sign of the null expansion does not depend 
on the scaling K: 9 > means expansion on average of the null generators, 
and 9 < means contraction on average. In Minkowski space, the future null cone 
^ = dl+{p) \ {p} has 9>0, and the past cone, ^ = dl~ip) \ {p}) has 9 < 0. 

It is useful to understand how the null expansion varies as one moves along a 
null generator of Let s — > 77(5) be a null geodesic generator of .yK, and assume 
K is scaled so that rj is affinely parameterized. Then it can be shown that the null 
expansion scalar 6 = 6{s) along 77 satisfies the propagation equation, 

(7.2) ^ = -Ric(,;',^')-^' ^ 



ds n — 1 

where cr > 0, the shear scalar, measures the deviation from perfect isotropic ex- 
pansion. Equation (7.2) is known in the relativity community as the Raychaud- 
huri equation (for a null geodesic congruence) [204], and, together with a timelike 
version, plays an important role in the proofs of the classical Hawking-Penrose sin- 
gularity theorems [204]. There are well-known Riemannian counterparts to this 
equation, going back to work of Calabi [76]. 

Equation (7.2) shows how the curvature of space-time infiuences the expansion 
of the null generators. We consider here a simple application of the Raychaudhuri 
equation. 

Proposition 7.1. Let ^ he a space-time which obeys the null energy condition, 
Kic{X,X) — RapX^X^ > for all null vectors X, and let JV he a smooth null 
hyper surface in ^ . If the null generators of jV are future geodesically complete 
then the null generators of -jy have nonnegative expansion, 9 >0. 

Proof. Suppose 9 < a.t p £ jV . Let 77 : [0,oo) — >■ jV , s — > rj{s), be the null 
geodesic generator of ^ passing through p = 7/(0); by rescaling K if necessary, 
we can assume rj is affinely parameterized. Let 9 = 9{s), s € [0,oo), be the null 
expansion of ^ along 77; hence 9{0) < 0. Raychaudhuri 's equation and the null 
energy condition imply that 9 = 9{s) obeys the inequality, 

(7.3) ? < - ^ - 



ds n — 1 

and hence 6* < for all s > 0. Dividing through by 9"^ then gives, 

which implies 1/9 ^ 0, i.e., 9 — > —00 in finite affine parameter time, contradicting 
the smoothness of 9. □ 
We wish to indicate the connection of Proposition 7.1 with the theory of black 
holes. In fact, this proposition is the most rudimentary form of Hawking's famous 
area theorem [204]. Let ^ be a standard black hole space-time, as defined for 
example in [204]. It is not necessary to go into the technical details of the definition. 
It suffices to say that in ^ there exists a region the black hole region, from 
which signals (future directed causal curves) cannot "escape to infinity" (recall the 
example of the Schwarzschild solution discussed in Section 3.1). The boundary of 
this region is the event horizon <S , which, in general, is a Lipschitz hypersurface 
ruled by future inextendible null geodesies, called its null generators. If S is smooth 
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and if its generators are future complete then Proposition 7.1 implies that S has 
nonnegative null expansion. This in turn implies that "cross-sections" of S are 
nondecreasing in area as one moves towards the future, as asserted by the area 
theorem. In the context of black hole thermodynamics, the area theorem is referred 
to as the second law of black mechanics, and provides a link between gravity and 
quantum physics. As it turns out, the area theorem remains valid without imposing 
any smoothness assumptions; for a recent study of the area theorem, which focuses 
on these issues of regularity, see [131]. 

7.2. Trapped and marginally trapped surfaces. We begin with some defini- 
tions. Let S = S"^^. n > 3, be a spacclike submanifold of co-dimension two in a 
space-time (^"+^, g). Regardless of the dimension of space-time, we shall refer to 
S as a surface, which it actually is in the 3 + 1 case. We arc primarily interested 
in the case where S is compact (without boundary) , and so we simply assume this 
from the outset. 

Each normal space of E, [TpE]^, p G S, is timelike and 2-dimensional, and 
hence admits two future directed null directions orthogonal to S. Thus, if the 
normal bundle is trivial, S admits two smooth nonvanishing future directed null 
normal vector fields /+ and /_, which are unique up to positive pointwise scaling, 
sec Figure 7.1. By convention, we refer to /+ as outward pointing and Z_ as inward 
pointing. In relativity it is standard to decompose the second fundamental form 




Figure 7.1. The nuU future normals l± to S. 

of S into two scalar valued null second forms x+ ^md x_, associated to 1+ and /_, 
respectively. For each p € S, x± ■ TpY, x TpS — > R is the bilinear form defined by, 

(7.5) X±{X,Y)=g{Vxl±,Y) for all X, y e T^S . 

A standard argument shows that x± is symmetric. Hence, x+ ^"^d X- can be traced 
with respect to the induced metric 7 on E to obtain the null mean curvatures (or 
null expansion scalars), 

(7.6) e± - tr^ x± = f'{x±\3 = divs^i . 

9± depends on the scaling of l± in a simple way. As follows from Equation (7.5), 
multiplying l± by a positive function / simply scales 6± by the same function. 
Thus, the sign of 6± docs not depend on the scaling of l±. Physically, (rcsp., 
0_) measures the divergence of the outgoing (rcsp., ingoing) light rays emanating 
from E. 

It is useful to note the connection between the null expansion scalars 9± and the 
expansion of the generators of a null hypersurface, as discussed in Section 7.1. Let 
jV^ be the null hypersurface, defined and smooth near S, generated by the null 

^-'^In many situations, there is a natural choice of "inward" and "outward". 
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geodesies passing through S with initial tangents /+ . Then 6+ is the nuU expansion 
of ^+ restricted to T,; 9- may be described similarly. 

Let M be a spacelike hypersurface in a space-time {^,g), with induced metric 
h and second fundamental form K , and suppose S is embedded as a 2-sided hyper- 
surface in A/. Then the null expansions 6± can be expressed in terms of the initial 
data (Af, /i, K) as follows. Since S is 2-sided in Af , it admits a smooth unit normal 
field v in M , unique up to sign. By convention, we refer to such a choice as outward 
pointing. Letting u denote the future directed unit normal to Af, = u + v (resp., 
L = u — I/) is a future directed outward (resp., future directed inward) pointing 
null normal vector field along E. Let 0± be the null expansion with respect to the 
null normal l±. Then, 

(7.7) 0± = ti-^K ± H , 

where tr-^K is the trace of the projection of K into E with respect to the induced 
metric on E, and H is the mean curvature of E in Af. 

For round spheres in Euclidean slices of Minkowski space, with the obvious choice 
of inside and outside, one has 6'_ < and 0+ > 0. In fact, this is the case in general 
for large "radial" spheres in asymptotically flat spacelike hypersurfaces. However, 
in regions of space-time where the gravitational field is strong, one may have both 
^_ < and 0-1- < 0, in which case E is called a trapped surface. For example 
the black hole region, < r < 2m, in [n + l)-dimensional Schwarzschild space- 
time (see Section 3.1) is foliated by spherically symmetric (n — l)-spheres, all of 
which are trapped surfaces. Under appropriate energy and causality conditions, 
the occurrence of a trapped surface signals the onset of gravitational collapse. This 
is the implication of the Penrose singularity theorem [317], the first of the famous 
singularity theorems. 

Theorem 7.2 (Penrose [317]). Let he a globally hyperbolic space-time with 
noncompact Cauchy surfaces satisfying the null energy condition. If contains a 
trapped surface E then ./# is future null geodesically incomplete. 

Recall from Section 7.1 that the null energy condition is the curvature require- 
ment, Ric (X, X) = Ra/sX^X^ > for all null vectors X. If a space-time ^ obeys 
the Einstein equation (2.8), then one can express the null energy condition in terms 
of the energy momentum tensor: ^ obeys the null energy condition if and only if 
TiX,X) = TijX'Xi > for aU nuU vectors X. 

In studying an isolated gravitating system, such as the gravitational collapse of a 
star, it is customary to model the situation by a space-time which is asymptotically 
flat. In this context, the assumption of Theorem 7.2 that the space-time admits 
a noncompact Cauchy surface is natural. The conclusion in the theorem of future 
null geodesic incompleteness is an indication that space-time "comes to an end" or 
develops a singularity somewhere in the causal future. However, the theorem gives 
no information about the nature of the singularity. 

Existence of vacuum asymptotically flat initial data sets, with one asymptotic 
region and containing compact trapped surfaces, has been established by Beig and 
6 Murchadha [38]. 

An intriguing question that arises in this context is whether a trapped surface 
can develop dynamically from initial data that did not contain any. This has been 
addressed by Christodoulou in [106], where a formidable analysis of the focusing 
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effect of sufficiently strong incoming gravitational waves is presented. Some further 
developments can be found in [250]. 

Remark on the proof of Theorem 7.2. Under the assumption of future null com- 
pleteness, one uses the trapped condition to show that the achronal boundary 
dI~^{Ti), which is generated by null geodesic segments issuing orthogonally from 
S, is compact. Then, flowing along the integral curves of a timclike vector field 
establishes a homemorphism from to any Cauchy surface S, contradicting 

the noncompactness assumption. □ 

To continue our discussion, consider again the general setting of a spacelike 
surface S""^ in a space-time with future directed null normal fields l± and 

associated null expansion scalars 6±. Focusing attention on just the outward null 
normal we say that S is an outer trapped surface (resp., weakly outer trapped 
surface) if 0-^, < (resp., 0+ < 0). If 0-)- vanishes, we say that E is a marginally 
outer trapped surface., or MOTS for short. In what follows we will be primarily 
concerned with properties of MOTSs. 

MOTSs arise naturally in a number of situations. As an outgrowth of their work 
on the positive energy theorem, Schoen and Yau [348] showed that suitable condi- 
tions on the energy density and momentum density of an asymptotically flat initial 
data set insure the presence of a MOTS; see also [379]. Next, as follows from our 
comments about the area theorem in Section 7.1, cross-sections of the event horizon 
in black hole space-times have nonnegative expansion ^ > 0. (By a cross-section, 
we mean a smooth intersection of the event horizon with a spacelikc hypersurface.) 
In the steady state limit this expansion goes to zero. Thus, it is a basic fact that 
cross-sections of the event horizon in stationary black hole space-times are MOTSs. 
For dynamical black hole space-times, MOTSs typically occur in the black hole re- 
gion, i.e., the region inside the event horizon. While there are heuristic arguments 
for the existence of MOTSs in this situation, based on looking at the boundary of 
the 'trapped region' [204, 369] within a given spacelike slice, using an approach put 
forth by Schoen [341], their existence has been rigorously established under natural 
physical conditions, flrst by Andersson and Metzger [10] in dimension three, and 
then by Eichmair [166] up to dimension seven; see Section 7.5. 

As noted earlier, MOTSs may be viewed as space-time analogues of minimal 
surfaces in Riemannian geometry. In fact, as follows from Equation (7.7), in the 
time-symmetric case {K = 0) a MOTS is simply a minimal surface in M. Of 
importance for certain applications is the fact, flrst discussed by Andersson, Mars 
and Simon [7, 8] , that MOTS admit a notion of stability analogous to that for 
minimal surfaces. 

7.3. Stability of MOTSs. In Riemannian geometry, a minimal surface (surface 
with vanishing mean curvature) is stable provided, for a suitable class of variations, 
the second variation of area is nonnegative, S^A > 0. Stability of minimal surfaces 
can also be characterized in terms of the associated stability operator. This latter 
approach extends to MOTSs, as we now describe. 

Let S be a MOTS in Af with outward unit normal v. Consider a normal variation 
of S in M, i.e. a map F : (-e, e) x E -J> M, such that (i) F(0,-) = ids and 
(ii) ^Ij^o = 01/, g C°°(E). Let 0{t) denote the nuh expansion of E* := F{t, E) 
with respect to It = u + vt, where u is the future directed timelike unit normal to 
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M and vt is the outer unit normal to in M. A computation shows 

80 

t=o 

where L : C°°(E) C°°{Y.) is the operator [8,190], 

(7.9) - + 2{X, \/cf>) + - (p + J(;.)) - + div X - |X| 

In the above. A, V and div are the Laplacian, gradient and divergence operator, 

respectively, on E, S* is the scalar curvature of S, p and J are the energy density 

and momentum density, respectively, as defined in Equations (5.3), (5.4), X is the 

vector field on E defined by taking the tangential part of V^u along E, and ( , ) 

denotes the induced metric on E. 

In the time-symmetric case, in (7.8) becomes the mean curvature H, the vector 

field X vanishes and L reduces to the classical stability operator (linearization 

of the mean curvature operator) of minimal surface theory. In analogy with the 

minimal surface case, we refer to L in (7.9) as the stability operator associated 

with variations in the null expansion 0. Although in general L is not self-adjoint, 

its principal eigenvalue (eigenvalue with smallest real part) Ai(L) is real. Moreover 

there exists an associated eigcnfunction which is positive on E. Continuing the 

analogy with the minimal surface case, we say that a MOTS is stable provided 

'^i(-^) ^ 0- (In the minimal surface case this is equivalent to the second variation 

of area being nonnegative.) It follows from basic properties of L that a MOTS E 

is stable if and only if there exists a normal variation of E, with > 0, such that 

^1 > 
at \t=o - 

Stable MOTSs arise naturally in physical situations, for example, as outermost 
MOTSs. We say E is an outermost MOTS in M provided there are no weakly 
outer trapped {0+ < 0) surfaces outside of, and homologous to, E. We say E is 
a weakly outermost MOTS in M provided there are no outer trapped {9+ < 0) 
surfaces outside of, and homologous to, E. Clearly, "outermost" implies "weakly 
outermost". Moreover we have the following [7,8]: 

Proposition 7.3. Weakly outermost MOTSs are stable. 

To see this consider a variation {Ej} of a weakly outermost MOTS E with 
variation vector field V ~ (pv, where is a positive eigcnfunction associated to the 
principal eigenvalue Ai = \i{L). If Ai < 0, then Equation (7.8) implies ff |j^q = 
Ai0 < 0. Since 0{Q) = 0, this implies that, for small t > 0, Et is outer trapped, 
contrary to E being weakly outermost. 

A standard fact in the theory of black holes is that, for black hole space-times 
obeying the null energy condition, there can be no weakly outer trapped surfaces 
contained in the domain of outer communications (the region outside of all black 
holes and white holes). It follows that compact cross-sections of the event horizon 
in stationary black hole space-times obeying the null energy condition are stable 
MOTSs. Moreover, results of Andersson and Metzger [9,10] provide natural criteria 
for the existence of outermost MOTSs in initial data sets containing trapped regions; 
see Section 7.5. 

Stable MOTSs share a number of properties in common with stable minimal sur- 
faces. This sometimes depends on the following fact. Consider the "symmetrized" 
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operator Lq : C°°(I]) ^ C°°(E), 

(7.10) Lo(</.) = -A0+ + JH)-i|xp)0. 

formally obtained by setting X = in (7.9). The key argument in [190] shows the 
following (see also [8,189]). 

Proposition 7.4. Ai(Lo) > Ai(L). 

In the next subsection we consider an application of stable MOTS to the topology 
of black holes. 

7.4. On the topology of black holes. A useful step in the proof of black hole 
uniqueness (see Section 3.9) is Hawking's theorem on the topology of black holes [204] 
which asserts that compact cross-sections of the event horizon in 3 + 1-dimensional, 
appropriately regular, asymptotically flat stationary black hole space-times obeying 
the dominant energy condition are topologically 2-spheres. As shown by Hawk- 
ing [203], this conclusion also holds for outermost MOTSs in space-times that are 
not necessarily stationary. The proof in both cases is variational in nature, and 
relies on the classical Gauss-Bonnet theorem. Developments in physics related to 
string theory have lead to an increased interest in the study of gravity, and in par- 
ticular black holes, in higher dimensions; see e.g. [170] for a recent review. The 
remarkable example of Emparan and Reall [169] of a 4 -I- 1 asymptotically flat sta- 
tionary vacuum black space-time with horizon topology S^xS"^, the so-called "black 
ring", shows that horizon topology need not be spherical in higher dimensions. This 
example naturally led to the question of what are the allowable horizon topologies 
in higher dimensional black hole space-times. This question was addressed in the 
papers of [189,190], resulting in a natural generalization of Hawking's topology 
theorem to higher dimensions, which we now discuss. 

Consider a space-time (^""'"^, .g), n > 3, satisfying the Einstein equations (not 
necessarily vacuum), and for simplicity assume that the cosmological constant van- 
ishes, A = 0. Let M be a spacelike hypersurface in which gives rise to the initial 
data set (M, h, K), as in Section 7.3. Recall from Equation (5.15) that the dominant 
energy condition holds with respect to this initial data set provided /c > | J| along 
M, where p and J are the energy density and momentum density, respectively, as 
deflned in Section 5.1 (but with A set to zero). 

The following result, obtained in [190], gives a natural extension of Hawking's 
black hole topology theorem to higher dimensions. Recall, a Riemannian manifold 
S is of positive Yamabe type if it admits a metric of positive scalar curvature. 

Theorem 7.5. Let be a stable MOTS in an initial data set {Ar\h,K), 

71 > 3. 

(1) If p > \ J\ along E then E is of positive Yamabe type. 

(2) If p > \J\ along E then E is of positive Yamabe type unless E is Ricci flat 
(flat if n = 2,3), X = ^'^'^ P + Ji'^) = along E. 

In the time-symmetric case. Theorem 7.5 reduces to the classical result of Schoen 
and Yau [345], critical to their study of manifolds of positive scalar curvature, that 
a compact stable minimal hypersurface in a Riemannian manifold of positive scalar 
curvature is of positive Yamabe type. 
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The key to the proof of Theorem 7.5 is Proposition 7.4, which, since S is assumed 
stable, imphes that Ai(io) > 0, where Lq is the operator given in (7.10). Now, 
in effect, the proof has been reduced to the Riemannian case. Consider S in the 
conformany related metric, 7 = (j)^-^"f, where 7 is the induced metric on S and cj) 
is a positive eigenfunction corresponding to Ai(Lo). The scalar curvatures S and 
S of the metrics 7 and 7, respectively, are related by (compare (5.11)) 



where for the second equation we have used (7.10). Since p + J(i^) > p — |J| > 0, 
we have that > 0. By further standard metric deformations, the scalar curvature 
of S can be made strictly positive, unless various quantities vanish identically. 

According to Theorem 7.5, apart from certain exceptional circumstances, a stable 
marginally outer trapped surface E in a space-time ^ obeying the dominant energy 
condition is of positive Yamabe type. Assume for the following discussion that S 
is orientable. Then, in the standard case: dim^ = 3 + 1 (and hence dimS = 2), 
Gauss-Bonnet tells us that if S is positive Yamabe then E is topologically a two- 
sphere, and we recover Hawking's theorem. In higher dimensions, much is now 
known about topological obstructions to the existence of metrics of positive scalar 
curvature. While the first major result along these lines is the famous theorem 
of Lichnerowicz [265] concerning the vanishing of the A genus, a key advance in 
our understanding was made in the late 1970s and early 1980s by Schoen and 
Yau [343,345], and Gromov and Lawson [195,196]. Let us focus on the case: 
dim ^ = 4-1-1, and hence dim E = 3. Then by results of Schoen- Yau and 
Gromov-Lawson, in light of the resolution of the Poincare conjecture, S must be 
diffeomorphic to a finite connected sum of spherical spaces (spaces which are covered 
by the 3-sphere) and x S^^s. Indeed, by the prime decomposition theorem, E 
can be expressed as a connected sum of spherical spaces, 5^ x S'^'s, and K{-k, 1) 
manifolds (manifolds whose universal covers are contractible). But as E admits 
a metric of positive scalar curvature, it cannot have any isr(7r, l)'s in its prime 
decomposition. Hence the basic horizon topologies in dim = 4+1 are 
and S"^ X 5^ (in the sense that E is "built up" from such spaces), both of which 
are realized by nontrivial black hole space-times. It remains an interesting open 
question which topologies of positive Yamabe type can be realized as outermost 
MOTSs; see [351] for examples involving products of spheres. 

A drawback of Theorem 7.5 is that it allows certain possibilities that one would 
like to rule out: for example, the theorem does not rule out the possibility of a 
vacuum black hole space-time with toroidal horizon topology. (This borderline 
case also arises in the proof of Hawking's theorem). In fact, one can construct 
examples of stable toroidal MOTSs in space-times obeying the dominant energy 
condition. Such MOTSs cannot, however, be outermost, as the following theorem 
asserts. 
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Theorem 7.6 ( [189]). Let S be an outermost MOTS in the spacelike hypersur- 
face M, and assume the dominant energy condition (3.15) holds in a space-time 
neighborhood ofT,.^^ Then E is of positive Yamabe type. 

As an immediate corollary, we have that compact cross- sections of event horizons 
in regular stationary black hole space-times obeying the dominant energy condition 
are of positive Yamabe type. In particular, there can be no toroidal horizons. 

Theorem 7.6 is an immediate consequence of the following rigidity result. 

Theorem 7.7 ( [189]). Let T, be a weakly outermost MOTS in the spacelike hy- 
persurface M , and assume the dominant energy condition holds in a space-time 
neighborhood ofT,. If T, is not of positive Yamabe type then there exists an outer 
neighborhood [/ « [0, e) x E o/ S m M such that each slice = {t} x T,, t £ [0, e) 
is a MOTS. 

Theorem 7.7 is proved in two stages. The first stage, and the main effort, is to 
establish Theorem 7.7 subject to the additional assumption that M has nonpositive 
mean curvature, r < 0. This is a purely initial data result, and the proof is carried 
out in two steps. The first step involves an inverse function theorem argument to 
show that an outer neighborhood of E can be foliated by surfaces Et of constant 
null expansion, 0{t) = Cf. This uses the stability of E in a critical way. In the 
second step it is shown that all of these constants are zero, cj = 0. It is here where 
the sign of the mean curvature of M needs to be controlled. Once having proved 
Theorem 7.7, subject to the condition, t < 0, the next stage, which is actually easy, 
is a deformation argument (specifically, a deformation of M near E) that reduces 
the problem to the case r < 0. 

We remark that Riemannian versions of Theorem 7.7 had previously been con- 
sidered in [74, 75]. 

7.5. Existence of MOTSs. As mentioned earlier, compact cross-sections of the 
event horizon in regular stationary black hole space-times are necessarily MOTSs. 
In dynamical black hole space-times, it is typical for trapped or outer trapped 
surfaces to form in the black hole region. But the occurrence of an outer trapped 
surface in a spacelike hypersurface that obeys a mild asymptotic fiatncss condition 
leads to the existence of a MOTS, as follows from the existence result alluded to 
near the end of Section 7.2, and which we now discuss [10, 166,341]. 

Let be a relatively compact domain in a spacelike hypersurface M, with smooth 
boundary dfl. We assume that dfl decomposes as a disjoint union of components, 
dfl ~ Einn U Eout, where we think of Ejnn as the "inner" boundary and Eout as the 
"outer" boundary of il in AI. We choose the normal along Ejnn that points into £7, 
and the normal along Eout that points out of fi, so that both normals point towards 
the region exterior to Eout- Thus, Einn is outer trapped ii 9 < with respect to the 
future directed null normal field along Ejnn that projects into f2. We say that Eout 
is outer untrapped ii 9 > with respect to the future directed null normal along 
Eout that projects out of fl. Heuristically, if Eout is lying in a region where M is 
"flattening out" then we expect Eout to be outer untrapped. 

With this notation and terminology we have the following existence result for 
MOTS. 

^■^Note that since wc arc assuming A = here, the dominant energy condition is equivalent 
to the condition, G pi/Xfy^ > for all future directed causal vectors X, Y , where Gp,^ is the 
Einstein tensor. 
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Theorem 7.8 ( [9, 166]). Let A/" be a spacelike hypersurface in a space-time 
, g) , with n < 7. Let Q, he a relatively compact domain in M, with smooth 
boundary dQ ~ UEout; such that the inner boundary Ei„„ is outer trapped and 
the outer boundary Eont is outer untrapped, as described above. Then there exists 
a smooth compact MOTS in homologous to S]i„„. 

Moreover, it has been shown that the MOTS constructed in Theorem 7.8 is 
stable [5,9]. 

In the time-symmetric case, Theorem 7.8 reduces to a well-known existence result 
for stable minimal surfaces. In the time-symmetric case the barrier conditions in 
Theorem 7.8 simply say that is a mean convex domain. One can then minimize 
area in the homology class of a surface in J7 parallel to Sinn and apply standard 
compactness and regularity results of geometric measure theory [174] to obtain a 
smooth (provided dimM < 7) stable minimal surface in fl homologous to Sinn- 
However, since MOTS do not arise as stationary points of some elliptic functional, 
such a procedure docs not work for general initial data. A completely different 
approach must be taken. 

In fact, the proof of Theorem 7.8 is based on Jang's equation [238] which is 
closely related to the MOTS condition 9^ = 0. Given an initial data set (Af , g, K), 
consider the Riemannian product manifold, M = R x M, g = dt^ -\- g, and extend 
i^T to M by taking it to be constant along the t-lines. Given a function / on M, 
consider its graph, Nf ~ graph / = {{t,x) £ M : t ~ f{x),x e M}, equipped with 
the induced metric. Then Jang's equation is the equation, 

(7.11) iH + P){f) = H{f) + P{f)^0, 

where H{f) = the mean curvature of Nf and P = tr n^K (compare with Equation 
(7.7)). 

Schoen and Yau [347] established existence and regularity for Jang's equation 
with respect to asymptotically flat initial sets as part of their approach to proving 
the positive mass theorem for general, nonmaximal, initial data sets. In the process 
they discovered an obstruction to global existence: Solutions to Jang's equation 
tend to blow-up in the presence of MOTS in the initial data {M,g,K). This 
problematic blow-up behavior that Schoen and Yau had to contend with has now 
been turned on its head to become a feature of Jang's equation: In order to establish 
the existence of MOTS, one induces blow-up of the Jang equation. This is the 
approach taken in [10, 166] in a somewhat different situation. 

In order to obtain solutions to Jang's equation one considers the regularized 
equation. 

(7.12) {H + P){ft) = tft . 

In [166] Eichmair uses a Perron method to obtain solutions ft to (7.12) for t suf- 
ficiently small, with values that tend to infinity in a small collared neighborhood 
of Eout, and that tend to minus infinity in a small collared neighborhood of Sinn, 
as f — >■ 0. The construction of Perron sub and super solutions makes use of the 
barrier conditions. Using the "almost minimizing" property [166] of the graphs 
Nt := graph /t, one is able to pass to a smooth subsequential limit manifold A^, 
bounded away from 9f2, each component of which is either a cylinder or a graph 
that asymptotes to a cylinder. The projection into A/ of such a cylinder produces 
the desired MOTS. 
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Under the barrier conditions of Theorem 7.8, Andersson and Metzger [10] were 
able to make use of this basic existence result, to establish for 3-dimensional initial 
data sets, the existence of an outermost MOTS, thereby providing a rigorous proof 
of a long held "folk belief in the theory of black holes that the boundary of the 
so-called outer trapped region in a time slice of space-time is a smooth MOTS. 
The key to proving the existence of an outermost MOTS is a compactness result 
for stable MOTS which follows from the extrinsic curvature estimates obtained 
by Andersson and Metzger in [9], together with an area bound obtained in [10]. 
This area bound follows from an outer injectivity radius estimate for a certain 
class of MOTS, established by an interesting surgery procedure. The powerful 
methods developed by Eichmair in [166, 167] can be used to extend this result to 
n-dimensional initial data sets, 3 < n < 7. 

Consider, now, a foliation {Mt}, a < t < b oi a. region of a space-time (^"^^, 
by spacclike hypersurfaces Mt- If each slice Mt admits an outermost MOTS Sj, 
then the family of MOTSs {St} may form a hypersurface in space-time. A smooth 
spacclike hypersurface iJ in a space-time foliated by MOTS is called a 

marginally outer trapped tube (MOTT). In [7,8], Andersson, Mars and Simon have 
obtained a rigorous existence result for MOTTs. Consider a spacelike foliation 
{Mt}, a < t < b, and suppose Sj^ is a MOTS in Mt^. They prove that if St,, is 
strictly stable, i.e. if the principal eigenvalue of the associated stability operator 
(7.9) is strictly positive, then there exists a MOTT H such that for t close to to, 
St := n Mt is a MOTS. For further related resuhs see [6]. 

Many challenging questions concerning the global existence and behavior of 
MOTTs remain open. Building on the analysis of Dafermos [152], Williams [376] 
gave sufficient conditions on spherically symmetric black hole space-times satisfying 
the dominant energy condition, insuring that a spherically symmetric MOTT ex- 
ists, is achronal, and is asymptotic to the event horizon. Understanding the generic 
asymptotic behavior of MOTTs in the nonspherically symmetric case remains an 
interesting open problem. 

A MOTT satisfying certain supplementary conditions gives rise to the notion 
of a dynamical horizon, which provides an alternative, quasi-local description of a 
black hole. For an extensive review of dynamical horizons and related concepts, 
including physical applications, see [17]. 



Appendix A. Open problems 

We compile here a list of interesting open problems discussed in the paper. They all 
appear to be difficult, of varying degrees of difficulty, with some most likely intractable in 
the foreseeable future. 

(1) Classify all vacuum near-horizon geometries with compact cross-sections, i.e., Rie- 
mannian metrics on compact manifolds satisfying (3.19), p. 19, in low dimensions. 

(2) Remove the hypotheses of analyticity, non-degeneracy, and connectedness in the 
black-hole uniqueness Theorem 3.4, p. 23. 

(3) Construct a five-dimensional, stationary, /^-regular vacuum black hole with two- 
dimensional group of isometries, or show that there are no such black holes. More 
generally, classify such black holes. Compare p. 23. 

(4) Prove that the positive energy theorem holds in all dimensions without the spin 
assumption. Compare p. 24 and p. 37. 
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(5) Find the optimal differentiability conditions for Theorem 4.1. More generally, 
construct a coherent local well posedness theory for the evolution and constraint 
equations for metrics with low differentiability; compare Remark 4.2, p. 28. 

(6) Describe in a constructive way the set of solutions of the vacuum constraint equa- 
tions on compact, asymptotically flat, and asymptotically hyperbolic manifolds, 
with arbitrary smooth initial data. Compare p. 35. 

(7) Find a well posed initial boundary value problem for the vacuum Einstein equa- 
tions which is well suited for numerical treatment, and prove numerical conver- 
gence. See Section 4.6, p. 32. 

(8) Formulate, and prove, a precise version of Conjecture 6.10, p. 58. More generally, 
formulate and prove a precise version of the BKL conjecture, or find an open set 
of metrics developing a singularity which do not exhibit a BKL-type behavior. 

(9) Show that generic Bianchi IX orbits have a dense tj-limit set on the Kasner circle; 
compare p. 56. 

(10) Show that uniqueness and existence of maximal globally hyperbolic developments. 
Theorem 6.2, p. 44, holds in an optimal weak differentiability class. This requires 
revisiting the whole causality theory for metrics of low differentiability. 

(11) Remove the condition of closed generators in Theorem 6.3, p. 45. More generally, 
show that existence of Cauchy horizons, not necessarily compact or analytic, 
implies existence of local isometrics, or construct a counterexample. 

(12) Generalize Ringstrom's Theorem 6.4, p. 46, to Gowdy models on S'^ x , on S'^, 
and on lens spaces L{p, q). 

(13) Show that some twisting U{1) x (7(1) symmetric vacuum models on T^' have 
mixmaster behavior. More generally, find an open set of such models with mix- 
master behavior. Even better, analyze exhaustively the asymptotic behavior of 
those models; see Section 6.2.2, p. 47. 

(14) Show that degenerate asymptotically flat spherically symmetric Einstein-Maxwell- 
scalar field solutions are non-generic; compare Section 6.2.3, p. 48. Analyze what 
happens in solutions in which trapped surfaces do not form. More generally, prove 
weak and strong cosmic censorship within this class of space-times. 

(15) Find an open set of U{1) symmetric metrics where the dynamics can be analyzed 
in the contracting direction. More generally, analyze exhaustively the dynamics 
of those models, see Section 6.4.1, p. 48. 

(16) Show that polyhomogeneous initial data in Theorem 6.7, p. 52, lead to solutions 
with polyhomogeneous behavior at null infinity. 

(17) Show that the Kerr solution is stable against small vacuum perturbations, com- 
pare Section 6.6, p. 53. 

(18) Determine what topologies of positive Yamabe type can be realized as outermost 
MOTSs in space-times satisfying the dominant energy condition, compare p. 66. 

(19) Prove initial data equivalents of Theorems 7.6 and 7.7, p. 67. 

(20) Describe the generic behavior near the event horizon of marginally trapped tubes 
in black hole space-times, see p. 69. 
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